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Abstract: We develop a general framework for computing the holographic 2-point functions 
and the corresponding conductivities in asymptotically locally AdS backgrounds with an electric 
charge density, a constant magentic field, and possibly non-trivial scalar profiles, for a broad class 
of Einstein-Maxwell-Axion-Dilaton theories, including certain Chern-Simons terms. Holographic 
renormalization is carried out for any theory in this class and the computation of the renormalized 
AC conductivities at zero spatial momentum is reduced to solving a single decoupled first order 
Riccati equation. Moreover, we develop a first order fake supergravity formulalism for dyonic renor¬ 
malization group flows in four dimensions, allowing us to construct analytically inhnite families of 
such backgrounds by specifying a superpotential at will. These RG flows interpolate between AdS 4 
in the UV and a hyperscaling violating Lifshitz geometry in the IR with exponents 1 < 2 ; < 3 
and 9 = z -\- 1. For 1 < 2 : < 2 the spectrum of fluctuations is gapped and discrete. Our hope 
and intention is that this analysis can serve as a manual for computing the holographic 1- and 
2-point functions and the corresponding transport coefficients in any dyonic background, both in 
the context of AdS/CMT and AdS/QCD. 
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1 Introduction and summary of results 

Holographic techniques have been applied in recent years to a wide range of quantum systems 
exhibiting strong dynamics, both in high energy and condensed matter physics. Even though the 
application of holography to most real-world systems is entirely phenomenological, holographic 
models have been very successful in capturing certain qualitative aspects of strongly coupled sys¬ 
tems. Indeed, such models provide an effective theory description for certain strongly coupled 
systems, completely on par with an effective field theory description, which, as is the case for 
non-BCS superconductors [1, 2], is often not known. 
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A class of phenomena that have attracted particular attention lately are phenomena occurring 
in the presence of a strong magnetic field, and considerable effort has been devoted to understanding 
such phenomena both within the field theory [3-12] and the holographic [13-27] frameworks. In 
particular, the focus in [7, 10, 11, 21] was on condensed matter phenomena, such as the Quantum 
Hall Effect [3-6, 15, 18, 19, 22, 28-32], while [9, 23, 27] study the effect of the magnetic field 
on heavy ion collisions and [8, 10, 13, 16, 17, 20, 26, 33] explore the Stark, chiral magnetic, chiral 
vortical, and Nernst effects. In fact, many of these phenomena have topological origin, which allows 
for a unified understanding of both condensed matter and heavy ion physics in strong magnetic 
fields [12]. 

Thermoelectric conductivities have been computed holographically from dyonic backgrounds 
at finite temperature in the pioneering works [34-36]. The magnetic field plays a crucial role in 
regulating the DC electric conductivity at finite charge density, removing the Drude peak that 
appears in any system with translation invariance. Translation invariance can also be broken in 
a number of other ways, leading to momentum relaxation (see e.g. [26] and references therein). 
Numerical dyonic backgrounds with running scalars at zero and finite temperature have also been 
considered recently in the study of the Fractional Quantum Hall Effect [22]. Finally, dyonic back¬ 
grounds with a mass gap provide the interesting possibility to study thermalization in heavy ion 
collisions in the presence of finite charge density and magnetic field, extending previous works 
studying thermalization in confining backgrounds as an initial value problem [37-40] . 

The purpose of the present work is twofold. Firstly, we put forward a general framework for 
systematically extracting holographic data from asymptotically locally AdS dyonic backgrounds in 
the presence of running scalars, based on a radial Hamiltonian formulation of the bulk dynamics. 
Although the general techniques we present here are not new, we feel that their power in sim¬ 
plifying the extraction of the holographic data has not been appreciated enough. In this respect 
we emphasize three aspects of the holographic dictionary: i) holographic renormalization and its 
importance to preserve the Ward identities, ii) n-point functions can be extracted directly from 
the bulk solutions, without evaluating the on-shell action, and iii) there are fundamental advan¬ 
tages in formulating the fluctuation equations as first order Riccati equations. All three aspects 
fundamentally rely on a Hamiltonian formulation of the dynamics. Indeed, there is a lot of millage 
one can get by exploiting the manifest symplectic structure of the Hamiltonian dynamics. This 
is because the local RG formulation of a quantum field theory [41], which is the natural language 
in the context of holography, inherently equips the space of local couplings and operators with a 
symplectic structure. 

As we shall see, correctly renormalizing the theory is crucial for holographically deriving Ward 
identities like (2.41) that relate 2- and 1-point functions, and which are required for expressing all 
thermoelectric conductivities in terms of the electric ones only. Moreover, the derivation of such 
identities becomes almost a triviality in the Hamiltonian language, since they follow directly from 
the first class constraints reflecting local symmetries in the bulk theory. Point ii) is important here 
as well. Namely, the 1-point functions are identified with the renormalized canonical momenta, in 
the presence of sources. This has a very concrete consequence which leads to drastic simplification 
of the computation of e.g. the 2-point functions: the 2-point functions can be computed by fur¬ 
ther differentiation of the 1-point functions, without evaluating the on-shell action. This is linear 
response theory applied directly to the bulk. Point iii) highlights the fact that the Riccati equa¬ 
tions are precisely equations for the response functions directly, eliminating the sources from the 
problem at the very beginning. The advantages of doing this are that a) it becomes almost trivial 
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to correctly renormalize individual 2-point functions, computing only the terms that contribute to 
a given correlator [42], b) the numerical solution of the fluctuation equations is considerably sim¬ 
pler since only a single regularity condition must be imposed in the IR and the arbitrary sources 
have already been eliminated from the problem [43], and c) the fact that the Riccati equations 
follow from the Hamiltonian formulation of the dynamics implies that the response functions are 
automatically compatible with the symplectic structure of the theory, thus avoiding the necessity 
of evaluating the symplectic form to correctly identify the modes in generic coupled systems [44], 
The second aim of this paper is to present exact analytic dyonic backgrounds with running 
scalars, some of which are confining. Most of the holographic backgrounds that have been used to 
study phenomena at finite charge density and magnetic field are known only numerically, or in some 
asymptotic limit. This is particularly the case for backgrounds that involve running scalars. Here, 
we present infinite classes of new analytic zero temperature solutions corresponding to dyonic 
renormalization group (RG) flows. In fact, utilizing the radial Hamiltonian formulation of the 
dynamics we derive a fake supergravity description of such backgrounds, both at zero and finite 
temperature. In the zero temperature case, and in the context of a bottom up model, this allows 
us to generate new solutions simply by specifying a superpotential at will. 


The model 

The theory we consider is a generic bottom-up Einstein-Maxwell-Axion-Dilaton (EMAD) model 
described by the action 

5 = ' 


d''+^xV=^ {R[g] - - Z{^)d^xd^x - V{cP, x) - + Sgh + ^05, 

( 1 . 1 ) 


Jx4 

where = 8ttG is Newton’s constant and 


1 


Sgh = 7^ / d'^x\/^ 2A, 

JdM 


( 1 . 2 ) 


is the Gibbons-Hawking term ensuring that this action admits a Hamiltonian description. Moreover, 
Sgs stands for a Chern-Simons term whose explicit form depends on the spacetime dimension. For 
d = 3, which we will mostly focus on, the Chern-Simons term takes the form^ 




(1.3) 


All four functions H((/>, x), 5]((/)), Z^c^) and n(x) specifying the action are a priori completely 
general, provided that the theory admits asymptotically AdS solutions. These functions can be 
restricted by demanding that the action (1.1) can be obtained as a consistent truncation of a gauged 
supergravity action, or by imposing a specific global symmetry, such as SL{2, Z) invariance as e.g. 
in [22]. The equations of motion following from (1.1) take the form 


df^xduX - ^9piydpxd'’x^ - 2S (^F^pF^ p - ^g^yFp^F^^'^ = 0, 


2Ugcl> -V^- Z^dpxd^x - ^4>FpaFP^ = 0 , 

2Vp {Zd^x) - V ^- li ^ eP '^ P ^ Fp . Fp ^ = 0, 

Vp i^FP’^ + UeP’^P^Fp^) = 0 . 


(1.4a) 

(1.4b) 

(1.4c) 

(1.4d) 


^We write with where n is the Lorentz frame index corresponding to the 

direction orthogonal to the boundary of dM in A4. 
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Summary of results 

In this paper we focus on homogeneous background solutions of these equations of motion with a 
finite electric charge density, a constant magnetic field, and running scalars. The presence of the 
electric charge density and constant magnetic field imply that these background solutions break 
Lorentz invariance, even at zero temperature. One of our main results is a completely general 
description of such backgrounds in terms of a fake superpotential (3.14) and first order BPS- 
like flow equations (3.15). Since Lorentz invariance is broken even at zero temperature, the fake 
supergravity description we develop here applies equally to both zero temperature RG flows and 
finite temperature black hole solutions. For solutions with zero charge density and magnetic field 
it reduces to the standard fake supergravity for Poincare domain walls at zero temperature, but 
it provides a non-trivial generalization to black hole solutions. It should be emphasized that the 
fake supergravity description we develop here follows from Hamilton-Jacobi (HJ) theory and it 
is quite different from the first order formulation in e.g. [45]. In particular, the superpotential 
(3.14) depends in general on both the scalars and the warp factor. Moreover, once a solution of 
the superpotential equation is given, the solutions of the equations of motion follow immediately 
from the first order equations (3.15). This is quite distinct from the formulation of [45], where 
the superpotential effectively facilitates a change of variables in order to write the second order 
equations as a system of coupled first order equations. 

In the context of a bottom up model like the one we study here, where the functions specifying 
the action are a priori arbitrary, the HJ origin of our fake supergravity description is particularly 
useful. In particular, one can very easily construct solutions of the superpotential equation (3.14) 
which immediately lead to exact solutions of the equations of motion. We present such an infi¬ 
nite class of exact RG flows in four dimensions specified by an arbitrary function Wo{4>) in the 
superpotential (3.18) in Section 3. These RG flows interpolate between AdS 4 in the UV and a 
hyperscaling violating geometry with exponents 1 < z < 3 and 6 = z + 1 m. the IR. By computing 
the Schrodinger potential of the fluctuations we show that the subclass of flows with 1 < z < 2 has 
a discrete spectrum of fluctuations, while the spectrum is continuous for 2 < 2 : < 3. A complete 
analysis of these exact backgrounds will presented elsewhere [46]. 

Besides the description of general homogeneous backgrounds with finite charge density and 
constant magnetic field, our aim in this paper is to provide a general framework for computing the 
renormalized 2-point functions of the dual stress tensor, current and scalar operators in any such 
background and for a generic model of the form (1.1). In this respect we describe the holographic 
renormalization of any such theory in terms of a counterterm potential U{4>, x) that is obtained by 
solving a superpotential equation in a Taylor series in (p and y, and we explicitly write the renor¬ 
malized 1-point functions in the presence of sources in terms of the renormalized radial canonical 
momenta in (2.37) [47]. This identification not only provides general expressions for the VEVs 
in any dyonic background in (3.32), but also allows one to extract the 2-point functions directly 
from the canonical momenta, without the need to evaluate the on-shell action. Moreover, using 
this identification of the renormalized 1-point functions in terms of the canonical momenta enables 
us to immediately identify the first class constraints in the Hamiltonian formalism, which reflect 
local symmetries in the bulk, with the Ward identities (2.38), (2.39), and (2.40), reflecting global 
symmetries in the dual field theory. Since these Ward identities hold for the 1-point functions 
in the presence of arbitrary sources, differentiating with respect to the stress tensor and current 
sources one obtains two Ward identities for the 2-point functions, which are given in (2.41). As we 
show in Section 5.2, these Ward identities together with the Kubo formulas allow one to express all 
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thermoelectric conductivities in terms of the electric conductivities, generalizing the corresponding 
result of [36] for the dyonic Reissner-Nordstrom black hole to more general backgrounds including 
running scalars. 

Another important aspect of our general framework for computing the renormalized 2-point 
functions is that we formulate the fluctuation equations in terms of first order Riccati equations 
instead of second order linear equations. The advantage of the Riccati equations is that they 
are differential equations for the response functions directly, with only one integration constant 
per equation, which can be fixed by imposing an IR condition. This drastically simplifies both 
the holographic renormalization of the 2-point functions [42] and especially the numerical solution 
of the fluctuation equations, since the arbitrary source is factored out of the problem from the 
onset [43]. For the finite frequency and zero spatial momentum fluctuations we consider here 
the fluctuation equations reduce to a single Riccati equation (4.15) (or equivalently (6.41)), from 
which all conductivities can be determined. This Riccati equation generically is not integrable, but 
some progress can be made by considering various limiting cases. Particularly interesting is the 
small frequency expansion which we determine up to (and including) in Section 6.4. The 

corresponding result we obtain for the conductivities is given in (6.34) and agrees completely with 
that obtained in [34] for the dyonic Reissner-Nordstrom black hole. Moreover, we hnd that using two 
different Fade approximants based on the small frequency expansion of the response functions (6.33) 
captures completely both the hydrodynamic and small magnetic field approximations discussed in 
[36], including the location of the poles closest to the origin of the complex frequency plane. 

Organization of the paper 

The rest of the paper is organized as follows. In Section 2 we present the radial Hamiltonian formu¬ 
lation of the dynamics described by the action (1.1) and we set up the algorithm for holographically 
renormalizing any such theory. This leads to a holographic derivation of the Ward identities, which 
we later use to relate various 2-point functions and conductivities. In Section 3 we discuss a general 
class of homogeneous backgrounds with finite charge density, constant magnetic field and running 
scalars and we develop a fake supergravity formalism that allows us to construct an infinite family 
of RG flows in four dimensions. Moreover, general expressions for the renormalized vacuum ex¬ 
pectation values (VEVs) corresponding to these backgrounds are derived. Section 4 contains the 
analysis of the fluctuation equations around the backgrounds described in Section 3 for certain 
class of time dependent fluctuations at zero spatial momentum. In particular, after decoupling 
the fluctuation equations following [22, 34, 36], we reduce the system to a single decoupled first 
order Riccati equation [42, 43]. In Section 5 we show how the renormalized 2-point functions can 
be extracted directly from the radial canonical momenta and the solution of the Riccati equation, 
showing that evaluating the on-shell action is a completely redundant step. Moreover, we discuss 
how the various conductivities are related to the 2-point functions of the current and stress tensor 
through the Kubo formulas, and derive the relations implied by the Ward identities. Section 6 
addresses the solution of the Riccati equation in various limiting cases. In particular, we derive 
the UV and IR asymptotic expansions and obtain general perturbative solutions for small and 
large frequency. All the machinery developed in the previous sections is applied to the dyonic 
Reissner-Nordstrom black hole in four dimensions in Section 7, where we solve numerically the 
Riccati equation and compare with the various approximate solutions discussed in Section 6. We 
end with some concluding remarks in Section 8, and three appendices where we provide a few 
more examples of exact RG flows in various dimensions, as well as the Gauss-Codazzi and general 
fluctuation equations following from the action (1.1). 
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2 Radial Hamiltonian formalism and the holographic dictionary 

The on-shell action is usually the starting point for any holographic computation since this is iden¬ 
tified with the generating function of the dual field theory [48, 49]. However, unless one wants 
to evaluate the free energy of a particular state, corresponding to the value of the (renormalized) 
on-shell action, evaluating the on-shell action is a redundant step since the radial canonical mo¬ 
menta are holographically identified with the one-point functions of the dual operators [47], which 
contain the same information. In particular, the vacuum expectation values (VEVs) of the dual 
operators correspond to the values of the (renormalized) radial momenta on a particular solution, 
while higher-point functions are obtained by successively differentiating the canonical momenta 
with respect to the sources. All n-point functions can therefore be extracted directly from the 
bulk solutions, without the need to evaluate the on-shell action. This is especially important in 
Lorentzian signature where it is often much harder to evaluate the on-shell action as a function of 
the sources than the canonical momenta. Moreover, the radial Hamiltonian formulation of the bulk 
dynamics allows one to most straightforwardly carry out the procedure of holographic renormaliza¬ 
tion, either fully non-linearly [47, 50, 51] or perturbatively to the desired order that contributes to 
the re-point functions one wants to renormalize [42]. Finally, the holographic Ward identities follow 
immediately from the radial Hamiltonian dynamics since they correspond to first class constraints. 

In order to formulate the bulk dynamics in a Hamiltonian language we begin by decomposing 
the bulk variables in components along and transversely to the radial coordinate, as in the standard 
ADM treatment of gravity [52], except that the Hamiltonian time now is the radial coordinate 
instead of real time. In particular the bulk metric is written in the form 

ds^ = {N'^ + NiN^)dr‘^ + 2Nidrdx^ -|- jijdx^dx^, ( 2 . 1 ) 

in terms of the lapse function N, the shift function W and induced metric 'jij, while the t/(l) gauge 

field is decomposed as 

A = adr + Aidx''. (2-2) 

In terms of these variables the bulk Ricci scalar becomes 

R[g] = R[-f] +K^- KijK^^ + {-2Kn'^ + 2rfV^,n ^), (2.3) 

where Kij is the extrinsic curvature given by 

~ (Tj ~ DiNj — DjNi ), (2.4) 

and re^ = is the unit normal vector to the constant r hypersurfaces. Here and in the 

following a ’ denotes a derivative with respect to the radial coordinate r and Di represents the 

covariant derivative with respect to induced metric 'jij. Inserting these expressions in the action 
(1.1) and using the identities 

/ 1 _ Af* \ 

= iVV=7, , (2.5) 

\-§7 + 
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we find that the action can be written in the form S = J drL and Scs = / drLcs with 


- dia - N^Fki){Aj - dja - N^Fij) 

- Z{^)f^dad,X - - Vi^, X)} , (2.6) 

and provided d = 3, 

Lcs = I A^^.J^AIl{x)FF^Ai - dia)F^k- (2.7) 

This Lagrangian is the basis of the radial Hamiltonian analysis. 

The canonical momenta following from this Lagrangian are 

= ( 2 , 8 a) 

= <"■***■> 

»x=^ = -; 3 ^ 2 (((.){i-iV 8 ix), ( 2 , 8 c) 

(n(x)e«‘r,l) , ( 2 , 8 d) 


while the canonical momenta conjugate to the fields N, Ni, and a vanish identically since the 
corresponding velocities do not appear in the Lagrangian. Note that the only contribution of the 
Chern-Simons term is the last term in the canonical momentum conjugate to the gauge field Ai. 
Clearly this term is present only for d = 3. Using these expressions for the canonical momenta we 
can evaluate the Legendre transform of the Lagrangian to obtain the radial Hamiltonian, namely 

H = y d'^x (^TT^^'yij + '^4:4> + Ti'xX + '^^Ai'^ ~ ^ ^ J 

where 

'H = - -^= (^2 (^XikXji - 

+ \F,-^i<P) + ^^/^U{x)€i^^Fk|^ ^TT* + ^y/^UixyP'^Fpq^ + ^Z-\(P)7tI + 

+ {-Rb] + Z{(t>)d,xd\ + + V{b x) + , (2.10a) 

W = - (^TTj + ^^U{x)eFFk^ + Tr^d\ + (2.10b) 

F = -DiTT\ (2.10c) 

Since the canonical momenta for the fields N, Ni and a vanish identically, the corresponding 
Hamilton’s equations imply the constraints 

n = W = F = {). ( 2 . 11 ) 

These are first class constraints, reflecting the diffeomorphism and gauge invariance of the bulk 
theory, and imply that the Hamiltonian is identically zero on-shell. These first class constraints are 
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the starting point for the analysis of the bulk dynamics and for the construction of the holographic 
dictionary. 

As is well known from HJ theory, the first class constraints (2.11) together with the condition 
that the canonical momenta be expressible as gradients of a functional 5[7, A, x, </>]; i-e. 


5S 

i 55 

5S 



6Afi 

'^X ~ ’ 

5x 



( 2 . 12 ) 


are completely equivalent to the full set of second order equations of motion. In particular, inserting 
these expressions for the canonical momenta in the constraints (2.11) leads to a set of functional 
differential equations for the function S[’j,A,x,(p]- These are the HJ equations for Hamilton’s 
principal function S['y,A,x,<i)], which is also identified with the on-shell action evaluated with a 
radial cut-off. The most general solution of the second order equations of motion is described 
by a complete integral of the HJ equations, i.e. a solution that contains as many integration 
‘constants’ (in this context functions of the transverse coordinates) as generalized coordinates. 
These integration constants can be thought of as the ‘initial momenta’ - the renormalized momenta 
denoted with a ^ in (2.36)- and as we shall see are identified holographically with the renormalized 
1-point functions. 

Given a complete integral of the HJ equations one can immediately write down hrst order 
BPS-like flow equations by identifying the expressions (2.8) and (2.12) for the canonical momenta. 
With the gauge choice iV = 1, iVj = 0 and a = 0, which we will adopt from now on, these flow 
equations become 


iij — 


x = - 
A = - 




liklji 


d-l 


6S 




lijlki 


5S 

hki' 


-T.-\(t>)Ii{x)ei^^F,k. 


(2.13a) 

(2.13b) 

(2.13c) 

(2.13d) 


Integrating these first order equations leads to another set of integration constants, which are 
identified holographically with the sources of the dual operators. The resulting space of solutions 
is parameterized by 2n integration constants, where n is the number of generalized coordinates, 
and hence spans the full space of solutions of the second oder equations of motion. This observa¬ 
tion is practically very useful since not only it drastically simplifies the procedure of holographic 
renormalization, but it also often leads to new exact background solutions. 


2.1 Holographic renormalization 

Holographic renormalization [47, 50, 51, 53-55] can be understood as the systematic procedure for 
determining the boundary terms required to render the bulk variational problem at infinity well 
posed. Once the variational problem is well posed the finiteness of the on-shell action is automatic 
[44, 56]. In fact, the boundary term required is always a specific solution of the radial HJ equation, 
not only for asymptotically locally AdS spaces [50] , but even more generally [44] . This should not be 
surprising given that a complete integral S of the HJ equations coincides with the on-shell action 
evaluated on the most general solutions of the equations of motion. In particular, inserting an 
asymptotic complete integral in the flow equations (2.13) and integrating them asymptotically is a 
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very efficient way to derive the Fefferman-Graham asymptotic expansions and their generalizations 
for non asymptotically AdS spaces. The key object therefore in renormalizing the theory and 
constructing the holographic dictionary is the general asymptotic solution of the HJ equation. 

If the bulk theory is renormalizable, all the ultraviolet (in the QFT sense) divergent terms in 
this general asymptotic solution of the HJ equations must be local in transverse derivatives, since 
they correspond to the local counterterms required to cancel the ultraviolet divergences of the on- 
shell action. The constraints = 0 and T" = 0 respectively impose invariance under transverse 
diffeomorphisms and U{1) gauge transformations. It follows that as long as we look for a local, 
gauge and transverse diffeomorphism invariant functional S['y,A,x,(p], the only equation that we 
have to solve is the Hamiltonian constraint Ti = 0. 

In fact we only need to solve the Hamiltonian constraint asymptotically, but in a covariant way, 
i.e. every term should be a function of the induced fields, without any explicit radial dependence.^ 
This can be achieved by formally expanding the functional S['y,A,x,(p] in a covariant expansion 
in eigenfunctions of a suitable functional differential operator 6. For asymptotically locally AdS 
spaces this can be the dilation operator [47], but the form of this operator depends on the various 
functions of the scalars (p and x that parameterize the Lagrangian, since these determine the leading 
asymptotic behavior of the induced fields. If we do not specify the explicit form of these functions, 
or even if we want to consider non asymptotically AdS backgrounds, then a covariant expansion 
still exists but in eigenfunctions of the operator [57] 

6^= [ (2.14) 

In the absence of a vector field this operator counts transverse derivatives, but this is not strictly 
the case when there is a Maxwell field in the theory. What (2.14) always counts is the number 
of inverse induced metrics, 7 ®-^, which coincides with the number of derivatives in the absence of 
a Maxwell field. When the solution of the HJ equations can be expanded both in eigenfunctions 
of the dilatation operator and of 6^, then the two expansions are simple rearrangements of each 
other. For example, different powers of the scalars would appear at different orders in the expansion 
according to the dilatation operator while they would all be part of the zero order solution of the 
expansion in eigenfunctions of 6.y. 

Before we proceed with the recursive solution of the HJ equation T-L = 0 for the function S, a 
few comments are in order regarding the role of the Chern-Simons term. Note that using ( 2 . 12 ) we 
can write 

+ ;|%/=7n(x)e'"'i"fc/ = ^^/=^^'(x)e®^■'=A,afcX + ^ (‘^ + )^ / . 

(2.15) 

It follows that if and only if n(x) is a constant then the effect of the Chern-Simons term is simply 
a shift of the on-shell action S according to 

S^S = S + ^j d3x7=^n(x)e®^'^AiF,fc, (2.16) 

where S satisfies the HJ equation derived from a Lagrangian without a Chern-Simons term. Cru¬ 
cially, this boundary term is gauge-invariant if and only if n(x) is a constant. In that and only 

^Keeping the boundary dimension d arbitrary allows one to write the conformal anomaly in covariant form as 
well, without any explicit dependence on the radial cut-off. However, the cut-off dependence emerges via dimensional 
regularization by replacing the pole in the coefficients of the asymptotic solution of the HJ equations with the radial 
cut-off [47]. 
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that case therefore we can simply determine S by solving the Hamiltonian constraint in the ab¬ 
sence of a Chern-Simons term and then replace S with S. For non-constant n(y), however, the 
Hamiltonian constraint in the presence of the Chern-Simons term cannot be reduced to that with¬ 
out such a term. Hence, to determine S one must solve directly the Hamiltonian constraint in 
(2.10). Nevertheless, we only need to determine the divergent part of the on-shell action S and 
it is straightforward to show that provided the equations of motion admit asymptotically locally 
AdS solutions the Chern-Simons term is asymptotically finite.^ For the purpose of determining the 
boundary counterterms therefore, even for non-constant n(x), we can always solve the Hamiltonian 
constraint corresponding to the theory without a Chern-Simons term. As we shall see though, the 
CS term leads to a gravitational anomaly in the holographic Ward identities. 

In order to construct the asymptotic solution of the HJ equation for general potentials in the 
Lagrangian we therefore proceed by formally expanding S as 


— 5(0) -I- 5(2) + • • • , (2-17) 

where each term 5(2fc) in this expansion is an eigenfunction of 5^ with eigenvalue d — 2k, i.e. 

^■y<S{2k) = {d — 2/i;)5(2fc)- (2-18) 

The zero order term 5(o) must contain no transverse derivatives and be gauge invariant and so it 
must take the form 

‘^(0) ^ ^ J (2-19) 

for some function of the scalars U{(p, x)- Inserting this in the Hamiltonian constraint % = Q leads 
to the ‘fake superpotential’ equation 


Ul + Z-\^)Ul - ^1/2 = H(0, x), 


( 2 . 20 ) 


where the subscripts denote partial derivatives with respect to the corresponding field. In asymp¬ 
totically locally AdS spacetimes this PDF can be solved by seeking a solution U (</>, x) in the form 
of a Taylor expansion, with the quadratic term determining the scaling dimension of the dual 
operator. In particular, for asymptotically locally AdS spacetimes the potential takes the form 


ml(t? + m\x^ + ••• , 

where the scalar masses must satisfy the Breitenlohner-Freedman bound [59] 


( 2 . 21 ) 


2 r2 2 r2 \ ^ 

m^L , m^L > - ( - 


and are related to the dimensions and of the dual operators as 

=-A^{d-A^), mlL^ =-A^{d-A^) 


( 2 . 22 ) 


(2.23) 


The possible solutions of (2.20) are of the form 


ui^,x) = 


d-l 




(2.24) 


^This would not necessarily be the case for a massive gauge field though. See e.g. [58]. 
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where and ^y. can each take two values, respectively or d — and A^^, or d — A^^,. However, 
only a solution of the form 


u{4>,x) — ——I- + 


(2.25) 


can be used as counterterms [42] and therefore one must seek a solution of this form. 

Once the relevant solution U{4>, x) is determined by solving ( 2 . 20 ), the higher order terms 5 ( 2 fc) 
are determined through the linear equations obtained by substituting the formal expansion (2.17) 
in the Hamiltonian constraint, namely^ 


2k^ 


1 


1 


^27r(o)*7r(2fc)i - ^^^yY7r(o)7r(2fc) + :^'n'4>(o)'^4,(2k) + ^^x(0)%(2fc)J - '^(2fc) 


where 


TT 


{2k) 


<^^( 2 fc) 


dSfOL’') 

vr(2fc) =-^> ^(2A:)x 




(2k) 


6S 


Sxij ’ 5Ai ’ 6x 

and the inhomogeneous source of this linear equation is given by 

1--'''L7JW^V 

J2 


TT 


{2k)<i> 


{2k) 


6(t) 


(2.26) 

(2.27) 


T^{2) (-^[ 7 ] + Z{4>)dixd\ + , 


9 fc—1 
K.^ 

'^{2k) = - 


T.FiiF'-^ - -j= j^27r(2)*7r(2)^ - j^t^{2) + 7r(2)*7r(2)i + 2^<i>(2)^ + ^ 

2 \ 

y~l ( 27r(2£)j7r(2fc-2£)i — j _ (21)'^{2k-2t) + ^{2i) {2k-2t)i 

i=l ^ 

+ 2'^<f>{2e)^4>{2k-2e) + ^^x( 2 £)^x( 2 fc- 2 £)^ ) k > 2. 


2 I 7-1.^ 2 

-Z 7rx(2) 


(2.28) 


Using the explicit form (2.19) of the leading order solution 5(o) these recursive equations become 


^ 1 

“ ^ {2fc) + Ufp7r^(^2k) + Z Uy^TTy^[2k) ='Tl{2k)- (2.29) 

This can be simplified further by writing 

S{2k) = J d'^x£(2fc)) (2.30) 

and using the identity 

^lij'^{2k)^ + ^Ai'n:(^2k) + <5<^7J'0(2fc) + ^XT^xi^k) = ^-^(2fc) + diV{2k) ^ (2.31) 

for some vector field V[ 2 k) ■ Applying this identity to the variations under 5^ and using the freedom 

to dehne T( 2 A:) up to a total derivative we arrive at the relation 

27r(2fc) — [d ~ 2fe)T(2fc), (2.32) 

which allows us to simplify the recursion relations to 

("m ^ J ■ (^) = ’^(2-=)' (2-33) 


■‘As discussed earlier, we use the Hamiltonian constraint without a Chern-Simons term since this term only 
contributes at finite order in the on-shell action. 
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These linear functional PDEs can be solved systematically following the procedure developed in 
[57]. Importantly, only the inhomogeneous solution, which is unique, contributes to the divergences. 
The homogeneous solutions are ultraviolet hnite. The counterterms are therefore defined as 


[d/2] 

Set •= ~ 'y ^ S^2k)- 

k=0 


(2.34) 


Carrying out this calculation keeping the dimension d as a parameter one finds that certain terms 
contain a pole at particular dimensions, e.g. l/{d — 4). Such terms lead via dimensional regular¬ 
ization to cut-off dependence according to the replacement rule 


1 

d-4 


(2.35) 


where To is the radial cut-off in the canonical radial coordinate r [47, 57]. The sum of all such terms 
is then identified with the holographic conformal anomaly [53]. 

Given the local counterterms Set the renormalized action (evaluated at a radial cut-off) is given 
by 

Sren := S + Set = j d'^X 4>Tt^ + X^x) > (2-36) 

where the quantities vf®, and are arbitrary functions that correspond to integration 
‘constants’ of the HJ equation. They correspond to the renormalized canonical momenta and can 
be identified with the renormalized 1-point functions of the dual operators through the relations^ 


(2.37) 


It should be emphasized that these are the 1-point functions in the presence of arbitrary sources and 
so any higher-point function can simply be obtained from these by further functional differentiation. 
These expressions are central to our subsequent analysis, and more generally to the holographic 
dictionary, since they allow one to extract the n-point functions of the dual operators directly from 
the radial canonical momenta, i.e. from the bulk solution of the equations of motion, without 
having to evaluate the on-shell action in terms of the sources. 


(rb) = 




13 


1 5Sr 


x/=7 6(j) 






TT 




(J®) = 


{Ox) = 


1 SSeen 
^/^ SAi 
1 dSep' 


1 


1 


TT 


\/=7 \/=7 


TTv 


2.2 Holographic Ward identities 

The identification of the renormalized 1-point functions with the renormalized canonical momenta 
in (2.37) allows us to translate the constraints 77® = 0 and T" = 0 into Ward identities for the 1-point 
functions in the presence of arbitrary sources. Since these constraints are linear in the canonical 
momenta, they hold at each order of the expansion (2.17) in eigenfunctions of the operator 5^. In 
particular, they hold for the renormalized momenta leading respectively to the diffeomorphism, 

D,{Ti) + Fij{,P) + {0^)diX + {0^)di(^ = -^U{x)e^’^^PjFki, (2.38) 

K- 

®To avoid cluttering the notation we do not differentiate between the renormalized 1-point functions evaluated 
at the cut-off, as e.g. in (2.37), and the limit obtained by multiplying these 1-point functions with the appropriate 
factor of the cut-off and sending the cut-off to infinity, as in e.g. (3.32). 
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and C/(l), 


A(J*) = 0, (2.39) 

Ward identities, which reflect symmetries that are always present in the bulk and so they hold 
irrespectively of the asymptotic form of the background. Notice that the Chern-Simons term intro¬ 
duces a gravitational anomaly in the dual theory, given by the RHS of (2.38). For asymptotically 
locally AdS backgrounds we also have the trace Ward identity (for a derivation see e.g. [47]) 

+ {d- + {d- A^){0^)d^cj> = A, (2.40) 

where A is the conformal anomaly, given by the coefficient of the cut-off dependent terms in Set- 
Since these Ward identities hold in the presence of arbitrary sources, we can differentiate them 
with respect to the sources to obtain constraints on higher-point functions. In particular, differen¬ 
tiating (2.38) with respect to the sources of the stress tensor and the current we get respectively 

I),(I]^'(x)r'='(x')) - 2Dj (jf (r')^(x))5W(x,x')) - (x))I?,<5('')(x,x') + (r"^(x))A5(")(x,x') 

+ F,,(x)(r'='(x') J^(x)) + = 

- 4n(x)e^>‘'F,,Fp,7^'<5('')(x,x') = 0, (2.41a) 

Kj 

D,{Ti{^)j\^)) + (J^(x)) - 5^5,) 5(‘')(x,x') + F,,(x)(J^(x)j'=(x')) + (0^(x)j"(x'))9a 

+ (0<^(x) J"(x'))5i</. = (dpi - dp,) + 2Fi,{pdpp) x'), (2.41b) 

where all derivatives are with respect to x, the covariant delta function is defined through 

<5(^)(x,x') = -^<5('^)(x - x'), (2.42) 

and we have used that 

(J^'(x)T^'(x')) = (r^'(x')J^'(x)) +7'=^5 (‘^)(x,x')(J^'(x)), (2.43) 

and similarly for the 2-point functions of the stress tensor and the scalar operators. We will revisit 
these two Ward identities for the 2-point functions in Section 5. 

3 Dyonic backgrounds 

Having established the holographic dictionary for the model (1.1), we are now interested in com¬ 
puting the holographic 2-point functions in backgrounds of the form 


ds% = dr'^ + [-f{r)df -h dx^ + dy'^) , 

(3.1a) 

Ab = a{r)dt -|- —{xdy — ydx), 

(3.1b) 

(j)B = XB = XB{r), 

(3.1c) 


where 77 is a constant background magnetic held, /(r) is the blackening factor and A{r) is the 
warp factor. The gauge held strength on such backgrounds is given by 

Fb = dAs = adr A dt + Hdx A dy. (3-2) 
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Inserting the ansatz (3.1) in the Gauss-Codazzi equations (B.2) we obtain the following set of 
equations for backgrounds of this form: 


1 


V 


+ 


2S 


A + AidA+-f-^fj + 


id + Zx^ = 0 , 

/ + / l^dA - ^/-V) - 4Se-2^ + 6?) = 0, 


^-2A 




-,-dA 


Vf 


dr i - 2nF^ = 0 . 


■^ + d^)=0, 

(3.3a) 


(3.3b) 


(3.3c) 

= 0, 

(3.3d) 

aH = 0, 

(3.3e) 


(3.3f) 


Note that due to the presence of the Chern-Simons term (1.3) these equations make sense either 
in any d provided II = 0, or for any 11 in d = 3. In what follows we keep 11 arbitrary and work in 
d = 3, but our results can be adapted to general d provided II is set to zero. 

In the following it will often be more convenient to work with an alternative radial coordinate 
defined through 

dr = (3.4) 

so that the background metric in (3.1) takes the form 

— f{u)dt‘^ + dx^ + dy^j . (3-5) 


ds% = 


du^ 

f{u) 


Denoting with a prime differentiation with respect to the radial coordinate u, the background 
equations for d = 3 become 


A" + 2A' (^A' + + a'2) = 0, (3.6a) 

2 [A" - A'2) + + Zx'l = 0, (3.6b) 

/" + 2A'f' - 4Ee-^^ + a'^) = 0, (3.6c) 

4’b + 2 (^A' + -/ (j)'^ — -f — -ZfjyXB ~ ^ = 0; (3.6d) 

^x'b + + -f x'b ~ 2 ^ + ^4>^'bX'b + 4n^e ^a'H = 0, (3.6e) 

(Sa' + 2ndd)' = 0. (3.6f) 

Notice that the last equation can be integrated in general leading to 

a' = ^iQ-2UH), (3.7) 


where the integration constant Q, as we shall see momentarily, determines the time component of 
the vacuum expectation value of the conserved 1/(1) current and so corresponds to the background 
electric charge density. 
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3.1 First order flow equations and the fake superpotential 


The radial Hamiltonian formulation of the dynamics we developed in Section 2 allows us to describe 
any solution of the second order equations (3.3) in terms of first order equations and a ‘fake 
superpotential’. In order to derive these first order equations we observe that for backgrounds of 
the form (3.1) the canonical momenta (2.8) (in the gauge N = 1, Ni = a = 0) become 


_ i)i (<5-. _ (1 + f)5ioSjo) + ^ 

71-0 = - 

K 

K 


(3.8a) 

(3.8b) 

(3.8c) 

(3.8d) 


Moreover, a variation of the on-shell action with respect to a metric of the form (3.1) takes the 
form 


5S = 7r“(57it -h TT°‘^6xab = ^^Att + f 


(3.9) 


and hence 

^ - l)r^A. (3.10b) 


Eliminating A from the first two equations in (3.3) and replacing the velocities with partial deriva¬ 
tives of 5(H,/,(/), X, a) using these expressions leads to the HJ equation for backgrounds of the 
form (3.1), namely 

(^) [^f^fi^A-dfSf)-Sl-Z-^Sl^+Veff{A,^,x) = 0, (3.11) 

where 

VeffiA, 0, x) = VicP, x) + + 2S-i(<^)e-2('’'-i)^ (g - 2i/n(x))" . (3.12) 

Notice that the background magnetic field and electric charge makes the effective potential depen¬ 
dent on the warp factor A in addition to the scalars (p and x- The dependence of 5 on / and a 
can be eliminated by the separable ansatz 

5 = -^ I d'^x (^e^^f/^W{A, cP, x) + 2Qa) , (3.13) 


where the term proportional to Q accounts for the canonical momentum conjugate to the time 
component a of the vector field and the fake superpotential W satishes the equation 


(3.14) 


Note that this ansatz is capable of providing an almost complete integral of the HJ equation, with 
“almost” referring to the fact that this ansatz does not contain an integration constant for the 
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generalized coordinate /. Inverting the momenta (3.8) and using the ansatz (3.13) for the on-shell 
action S leads to the advertised first order equations 



(3.15a) 

h-id-.r^’ 

(3.15b) 

i 

II 

■-s. 

(3.15c) 


(3.15d) 

d = _s-ie-(‘^-2)^/V2 (^Q _ 2HU{x)^ . 

(3.15e) 


As is guaranteed by the HJ construction (and can be checked explicitly using the background 
equations (3.3)), given a fake superpotential that satisfies (3.14) any solution of the first order 
equations (3.15) leads to a solution of the second order equations of motion (3.3). 

3.2 Exact families of dyonic backgrounds 

In fact, the superpotential equation (3.14) together with the first order equations (3.15) amounts to 
a (fake supergravity [60]) solution generating technique for dyonic backgrounds with running scalars, 
which becomes particularly powerful in the context of a bottom up model where the potentials 
defining the action (1.1) are a priori unspecified. To support this claim we provide two explicit 
examples here: a superpotential that gives the general dyonic Reissner-Nordstrom black hole in 
arbitrary dimension, and an infinite class of RG flows between AdS4 in the UV and a hyperscaling 
violating Lifshitz geometry in the IR. Additional solutions of the superpotential equation (3.14), 
including in other dimensions, are presented in Appendix A. 


Generalized dyonic Reissner-Nordstrom black hole 

The dyonic Reissner-Nordstrom black hole in arbitrary dimension can be obtained immediately 
from the flow equations (3.15), which imply that a universal consistent solution, independently of 
the form of the potentials as long as they admit asymptotically AdS solutions, is given hy cj) = x = 

In that case the superpotential equation (3.14) can be integrated to obtain 


kkRAr(A) 


d-l I 2L2e- 1(Q - 2nR)2e-2('^-iM 

{d-i){d-2) (d- i)(d-4) 


(3.16) 


where C is an integration constant and 11 is understood to be zero unless d = 3. Introducing the 
new radial coordinate u through® e~^ = u/L one can integrate the first order equations (3.15) to 
obtain the general form of the backgrounds corresponding to the superpotential (3.16), namely 


= ^ (/ ^{u)dv? — f{u)dt‘^ -|- dx^) , 
2L-2('^-2)s-i(0)(Q-2n(0)R)2 
' + - {d-l)(d-2) -" 

r-{d-3) _ 

a = ao + ^S-i(0)(g - 2U{0)H)u'^-\ 


2(d-l) 


2L-2s(0)R2 ^ 

{d-l){d-4f 




(3.17a) 

(3.17b) 

(3.17c) 


®This is a priori a different radial coordinate than the one introduced in (3.4), but for the Reissner-Nordstrom 
black hole the two coincide. 
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From these expressions we identify the integration constant C in (3.16) with the black hole mass, 
namely M = L~'^C, while the integration constant ao in the gauge field is identified with the 
chemical potential 


Exact dyonic RG flows 

For d = 3 (AdS4), the superpotential equation (3.14) admits a solution of the form 

W{A,4>) = Wo{cl))^l + q^e-^A, 
where the parameter q is defined through 

+ (Q - 2noR)2So 1 = q^L-^, 


(3.18) 


(3.19) 


and the scalar potentials V{(f)) and S(())) are determined in terms of a single arbitrary function 
Waicti) as 

V{4>) = < - 

+ {q- 2 UoHfj:-\cp) = ^ (w'^ + . (3.20) 

In order to solve explicitly the second equation for S((/)) we need to distinguish two cases, depending 
on whether the ratio \Q — 2IiQH\/\H\ is smaller or greater than 1. To do this it is convenient to 
introduce two parameters, Tg and through the identifications 

Q - 2noR ^ f coth(re/2), \Q - 2IioH\ > |R|, 

SoR ■ |tanh(r^/2), |Q-2noR| < |R|, ^ 

so that the corresponding expressions for S(ij!)) are 


|Q-2noR| > \H\ : 


1 


1 , 


1 , 


X hFo + xhF. ]+\ + -1T2 - L-4tanh" r, 


|Q-2noR| < \H\ : 

S(<^) = ^l2So( 1 + tanh2(r^/2)) (w'^ + + 

(3.22) 

Note that since IFo ~ —2/L + 0{(l)‘^) in the UV, the quantities under the square root are always non¬ 
negative. Moreover, since depends explicitly on the parameters Tg or Tm, these parameters (i.e. 
the ratio of the electric charge density to the magnetic field) must be considered as a specification 
of the theory - not of the solutions. 

For either of the two cases, inserting the superpotential (3.18) in the first order equations (3.15) 
allows us to obtain the explicit form of the metric in the form 
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where the warp factor is expressed in terms of the scalar via 


A = - 


1 

2 



(3.24) 


This geometry is asymptotically AdS 4 in the UV. Assuming that in the IR^ </>—)• oo and the function 
Wo{4>) behaves as 

lTo((/>) ~ A > 0, (3.25) 

the flow equations imply 

(3.26) 

Introducing the radial coordinate 

V = e-A+V2A)9i^ (3 27) 

in the IR the metric asymptotes to the hyperscaling violating Lifshitz geometry [61-68] 

with Lifshitz and hyperscaling violating exponents 


_ + 3/2 

“ A2 + 1/2’ 


6 = z + 1, l<z<3. 


(3.29) 


These are type Illb hyperscaling violating backgrounds for 1 < z < 2 and type IVb for 2 < z < 3 
according to the classification of solutions of the null energy conditions in [69]. It is interesting to 
note that for the purely electric solution [H = 0), these IR geometries correspond to solutions of 
an Einstein-Maxwell-Dilaton theory with a vanishing scalar potential [61, 68-70]. In the present 
case, however, the scalar potential is not zero! The reason for this apparent puzzle is that the 
scalar potential is asymptotically subleading in the IR in the purely electric case, and there are 
necessarily subleading terms, contrary to the case of pure exponential potentials. Another point 
we must emphasize is that, contrary to what is often claimed in the literature, these geometries 
are not singular, since the concept of a curvature singularity in the presence of a diverging scalar is 
inherently ambiguous. Indeed, as it was shown in [69], these geometries are perfectly well behaved 
in the ‘dual frame’ where they become Lifshitz and the singularity is completely absorbed in the 
diverging scalar. It follows that there is no need for applying Gubser’s criterion for this kind or 
IR geometries, in the same way that this criterion is not applicable to RG flows between two AdS 
fixed points. 

In the far IR, the radial coordinate v is related to the conformal coordinate u introduced in 
(3.4) as 


u ~ < 


6 2A^-1, 

p-c^u 


(u* — u) - 


U —)■ +00 A^ < |, 

u —>■ + 00 , A^ = 2 , 
U U-, A^ > 


(3.30) 


where cq is a non vanishing constant that depends on the charges. Gomputing the Schrodinger 
potential for the fluctuation equations (4.12) one finds that for both H ^ 0 and H = 0 cases it 
behaves in the IR as 


Vsch 


-h'l + -/14 — ho 


_ 2(2A^-1) 
2a2 + i , 


(3.31) 


^The case 4> —>■ —00 is completely analogous, with A < 0. 
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where ~ is non-zero. We therefore expect that the spectrum of fluctuations is discrete and 
gapped for > 1/2 (1 < 2 ; < 2), continuous and gapped for = 1/2 {z = 2), and continuous and 
ungapped for A^ < 1/2 (2 < 2 : < 3). A full analysis of these RG flows, including the possibility of 
embedding them in gauge supergravity will appear elsewhere [46]. 

3.3 Vacuum expectation values 

Using the identification (2.37) of the 1-point functions with the renormalized canonical momenta we 
can now evaluate the vacuum expectation values (VEVs) of the dual operators in these backgrounds. 
Combining (2.37) and (3.8) we find that the general form of the renormalized 1-point functions 
evaluated in backgrounds of the form (3.1) take the from 

{Ti,) = - hm (^((d-l)A {dij - (1 + mo6jo) + 

r—>-oo \ K \ Z J 

(OA = lim + , 

r—>-oo \ ^ ) 

{OA) = lim f —L^(^)^ + e-3Ay.-i/2^ctA ^ 

r^oo y ^ y 

(/) = 

To simplify these expressions we need the form of the counterterms Set and the corresponding 
canonical momenta tt/I, and evaluated In backgrounds of the form (3.1). These counterterms 
can be computed in full generality using the algorithm presented in Section 2.1, but this analysis 
can be simplihed by focusing on background solutions only. In particular, in Section 2.1 terms 
involving the held strength Fij were counted as derivative terms and do not appear in the leading 
order solution (2.19) parameterized by the function Since the held strength is constant in 

backgrounds of the form (3.1) however, we can take its effect into account by a modihed function 
U{A,4>,x) that satishes the same equation (3.14) as the fake superpotential W{A, (j),x), instead 
of the simpler equation (2.20). In fact, since we need to determine only the divergent part of 
U{A,(p,x) we can drop the term proportional to Q — 2HII in the effective potential (3.12) since 
this term would only affect U at order which is always subleading relative to The 

counterterm function U{A, cf), x) for general backgrounds of the form (3.1) therefore can always be 
obtained as a Taylor expansion solution of the following PDE: 

^ (d + dA) = V{<p, x) + (3.33) 

As we emphasized in Section 2.1, the relevant solution must have a Taylor expansion of the form 
(2.25). Erom the form of (3.33) we immediately deduce that the counterterm function [/((^, y) 
will be a function only of the scalars cj) and x as long d < 4 since in that case the term involving 
the background magnetic field contributes at subleading order relative to It follows that 

for d < 4 the function U{4>,x) can be obtained by solving the simpler equation (2.20). Eor the 
marginal case d = 4 this term will contribute to U with a coefficient that has a pole of the form 
l/(d — 4) (see (3.16)) which via the dimensional regularization described in Section 2.1 leads to 
a cut-off dependent divergence and a related conformal anomaly. Eor d > 4 the magnetic field in 
(3.33) will contribute to the ultraviolet divergences in the standard way. 


koko)j + 


2e 


- 3 A 


ct 


Vf 


(3.32a) 

(3.32b) 

(3.32c) 

(3.32d) 
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Finally, having determined the counterterm function U{(p, x) for backgrounds of the form (3.1) 
through (3.33), we obtain the renormalized 1-point functions 

{Tij) = mn e*/'- (((d - 1)A - [/) (% - (1 + /)«»«,») + (^ “ ^ j («« - 

= i _^lim e*/‘ ((»' + U) {Stj - (1 + f)S„Sja] + ^ {Wa + Va) (Sij - iffldjo)) , (3.34a) 


(0«) = -;(2 (ZX + C;,) = -i ^lim (ir, + C/,), (3.34c) 

= (3.34d) 

As was mentioned earlier, we see from these expressions that the integration constant Q corresponds 
to the conserved electric charge density of the background. Moreover, the energy density e is given 
by 

e := (T/) = -{Tu) = "4 ((^ " 1)^ " 4 , (3-35) 

r^oo \ / 

while the pressure density corresponds to (no index summation) 

V := (T,,) = -4 bm ((d - l)A - U + ^- . (3.36) 

r-)-oo \ 2j d — 1 j 


4 Fluctuation equations at zero spatial momentum 

We next consider linear fluctuations around general backgrounds of the form (3.1), and for arbitrary 
potentials V{(j), x), S((/)), Z((f>) and n(x). In order to be able to decouple the fluctuation equations 
(which we present in full generality for d = 3 in Appendix C) we follow [22, 34, 36] and consider 
only fluctuations that are independent of the spatial transverse coordinates and we set certain 
components of the fluctuations consistently to zero. 

In particular, denoting the most general fluctuations that preserve the gauge N = 1, Ni = a = 0 
by 

lij = iBij + hij, Ai = ABi + ai, (l) = 4>B + tp, X = XB + T, (4.1) 

with Sj = we switch off the fluctuations Sj = = Sy = Sx = (f = r = at = 0 and 

only keep the components ax = ax{r,t), ay = ay{r,t), Sf = Sf{r,t) and = S^{r,t). For such 
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fluctuations the only non-trivial equations are 

Einstein xt: 

(d‘^ + dr - Sf = -4Se-2^ (dd^ + He-^^dtUy) , (4.2a) 

(d^r + (si + ^/“V) dr + 4Se-2^/-id2^ 5* = 4Se-2^/-i (dd, + He-^^dtay) , (4.2b) 

Einstein yt: 

(52 + (si - - 4Se-^^i?2^ 5f = -4Se-2^ (dd^, - He-^^dta^) , (4.2c) 

(52 + (si + ^/“V) 5, + 4Se-2^/-id2^ 5* = 4Se-2^/-i (dd^ - He-^^dta,) , (4.2d) 

Einstein rx: 

= -4Se“2^ (i/d5]^ + HfcLy + adta ^), (4.2e) 

Einstein ry: 

dtSl = -4Se-2^ {-HaSf - Hfd^ + adtUy), (4.2f) 

Maxwell x: 

dr (s/-i/2e^ {dSf + fa,]) = Y.f-^l'^e-\d‘la^ + HdtSl) + 2Ii^XBdtay, (4.2g) 

Maxwell y: 

dr (s/-V2e^ (^aSf + fdy)') = Y.r^l'^e-\dlay - HdtSf) - 2U^Ma.. (4.2h) 

Note that (4.2a) and (4.2b) as well as (4.2c) and (4.2d) are trivially related since 

Si = -r^sf, si = -r^si (4.S) 

The remaining equations can be decoupled by introducing the complexified variables 

sf = Sf ± iSf, ht± = htx ± ihty, a± = ax± iuy, (4.4) 

so that we can write 

(df + (si - ^/“ V) sf = -4Se-2^ (dd± ^ iHe-^^dta± - H^e-^^Sf) , (4.5a) 

dtSf = 4Se“^"^ [±iH{dSf + fd±) — ddta±) , (4.5b) 

dr (^aSf + /d±)) = J:f-^/^e-^{dfa± ^ iHdtSf) ^ 2iIi^XBdta^, (4.5c) 

or after Fourier transforming in time {dt —)■ ioj) 

dr (e^^r^^^Sf^ = - 4 Se ^/-^/2 ^ ojHe-^^a± - H‘^e-^^Sf) , (4.6a) 

uSf = 4Se“^"^ {^H{dSf + fd±) — ujda±) , (4.6b) 

dr (^S/-i/ 2 gA (^^g± ^ ^ ± ojBSf) ± 2ujU^XBa±. (4.6c) 

Note that multiplying (4.6b) with taking the radial derivative and substituting (4.6c) in 

the resulting expression gives back (4.6a), which is therefore not independent, unless cu = 0. 
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In order to decouple these equations for w 7 ^ 0 we replace a± with the linear combinations 


= u)a± T HSf, 


(4.7) 


in terms of which (4.6) become 

udr = - 4 Se ^/-^/2 ^ ^ , (4.8a) 

(a;2 - 4Se-2^/772) 5^^ = 4Ee-^^ (±77/4 - wdf±) , (4.8b) 

^ ± 774)) = ^ 2ujII^xb£±, (4.8c) 

where we have used the last equation in (3.3) to obtain (4.8c). Substituting now the expression for 
4 from (4.8b) into (4.8a) leads to the two decoupled equations for £± 


4 ± 51 (w, 77)4 ±5 o(^^,±77)£1± = 0, 


where 


giiuj,H) =dr\og 


-1 


with 


fl'o(<^,±77) =/ ^^(Q - 4S4) =F 4(9^ log Se^/ ^ 

M 


£L = uj‘^ - 4S774e"^^. 


In terms of the u coordinate these equations read 


4 + hi{uj, 77)4 + ho{oo, ±H)£± = 0, 


where 


(4.9) 

(4.10a) 

(4.10b) 

(4.11) 

(4.12) 


/ii(a;, 77) log 4 / 17-4 (4.13a) 

ho{uj,±H) =/-2(17-4S/e-2V2)^4/-ia'5„log|Sa'l7-4 (4.13b) 

11 

Riccati form of the fluctuation equations 

Finally, these linear second order fluctuation equations can be expressed in first order Riccati form 
by introducing the response functions TZ± [42, 43], namely 


4 = 'Tl±£±, 


(4.14) 


so that (4.9) take the Riccati form 


(4.15) 

Since these are first order differential equations there is only one integration constant for each, 
which is fixed by imposing suitable boundary conditions in the interior of the geometry. This is 
one of the advantages of the Riccati formulation of the fluctuation equations, since they compute 
directly the response functions, without any dependence on the arbitrary sources. 


TZ± + TZ\ + 5i(a;, H)TZ± + go{^^, ±77) — 0. 
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5 Renormalized Green’s fnnctions, Kubo formnlas, and transport coefficients 

Another advantage of the Riccati form (4.15) of the fluctuation equations is that, as we now 
show, the 2-point functions can be determined in general in terms of the response functions TZ±. 
Therefore, solving the Riccati equations directly determines the 2-point functions, without the need 
to evaluate the on-shell action and then take functional derivatives. The Riccati formulation of the 
fluctuation equations implements linear response theory in the bulk. 

5.1 Holographic Green’s functions and Ward identities 

In order to unravel the relation between the response functions TZ± and the 2-point functions we 
need to determine the canonical momenta (2.8), which as we have seen in Section 2 are identified 
with the 1-point functions, to linear order in the fluctuations. For the fluctuations we considered 
in Section 4 the only non-trivial components of the canonical momenta to linear order in the 
fluctuations are 

+ 2is« + T [s'> - Shg) j , (5.1a) 

TT* = -^\/-7b 5] (dS** - {5l.ay - Sya^) , (5.1b) 

K/ \ / K 

whose only non-zero components are 

(1) 1 . . 

^h=- [Sf + ^ASf) , (5.2a) 

(1) 9 / ^ 

=^\/=^7b (s {aSf + /d±) T , (5.2b) 

where = 7r^±z7r^. Using the dehning relations (4.14) for the response functions 

TZ± we can express the velocities Sf and a± in terms of the fluctuations £± as 

Sf = 4Se“2^0“^ {±Hfn± - ua) £±, (5.3a) 

a± = n-^ (w7^± T 4RSe"2^d) £±. (5.3b) 

Hence, the only non-zero components of the momenta to linear order in the fluctuations take the 
form 

(1) 1 (1) 

'^h = {C^h^t± + G^a«±) 1 = ^J^lB {C^a^t± + Ctaa±) , (5.4) 

where 

C±^(r, uj) = ^ 756 - 2 ^ (±R/7^± - cad) - i) , (5.5a) 

C'ha(u {±Hfn± - cad), (5.5b) 

C'^a(uw) = -Uja) - , (5.5c) 

and we have used the fact that — 2IIH = —Q is a constant (see (3.7) and (3.4)). 

In order to renormalize these expressions for the canonical momenta linear in the fluctuations 
we must take into account the contribution of the boundary counterterms, which for d = 3 take 
the form 

<Sct = -^ j d^x 7=7 ([/(<)), x)©((/), x)R[7]) , (5.6) 
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where U is determined through (2.20), while 0 is determined by the algorithm described in Section 
2.1. However, for backgrounds of the form (3.1), where the Ricci scalar of the induced metric 
vanishes, and for fluctuations of the form discussed in Section 4 it is straighforward to check that 
the counterterm involving the Ricci scalar does not contribute (see last equation in (C.2)). In 
particular, when expanded to second order in the fluctuations the counterterms for d = 3 take the 
form 

(2) 1 c r 

S'd = ^ J d^x J du (5.7) 

These counterterms lead to the renormalized response functions 


r\j 

= {±Hfn± -ujd)-A + , (5.8a) 

ClT^ir,uj) = (r,u;) = (±77/77± - cud), (5.8b) 

C:T^{r,uj) = CtAr.u:) = (sQ-^o; (±77/77± - cud) - . (5.8c) 


Notice that only the coefficients get renormalized. 

Using (2.37) we can now express the renormalized 1-point functions in terms of the renormalized 
momenta, namely 


(T*^) = (r*^) ± i{T^y) = —, 

yJ-lB 

(j±) = (J") ± i{jy) = (5.9) 

V-7i? 


or, after lowering the indices 


{Tt±) = 


V-lB 


„4:A jl/2;^± 


= 2e" 


h ’ 


{J±) = 


V-lB 


e 


= e-^f 




(5.10) 


The 2-point functions can now be obtained by functional differentiation of these 1-point functions 
with respect to the fluctuations. The precise coefficients of the functional derivatives with respect 
to ht± and a± that correspond to an insertion of respectively and Jzp can be determined by 
computing the variation 

SSren = 2tt^^ 6hxt+2TTy^6hyt+^^duy = ^ [2^'^6ht- +2TfJ^5ht+ + Td^5a- + Tf~5a+) , (5.11) 

where the factor of 2 in the metric momenta is due to the fact that we must sum over the two 
possible index combinations. It follows that 

(T,±{a,)) = (5,12) 

dht^[—uj) dazf(—uj) 

We now have all the ingredients in order to evaluate the renormalized 2-point functions. Com¬ 
bining (5.9) and (5.4) we get 


{r,±)= (ciri>,± + c^rat), {j±)= (c^rh,i +ci"%±) , 


(5.13) 














From (5.12) then follows that the 2-point functions are obtained as 


ont+[io) 

(r,+ (a,)J_(-a,)) = 2e-"‘j^(r,+ (L,)> = 2e^^C+:‘’‘(w), 

(J+MJ_(-a,)) = 2e-'‘;p^(J+M) = 2e"'cr“M, (5,14) 

da+fcjj 

and similarly for the fluctuations ht- and a_. In particular, the full set of renormalized 2-point 
functions that can be computed with the fluctuations we considered are 

{Tt+{u)Tt-{-u;)) = 2^1im a;)) , (5.15a) 

{TUu^m+i-co)) = 2^1im , (5.15b) 

(r*+(a;)J_(-a;)) = 2 lim (e5^/^C+;-(r,cu)) , (5.15c) 

(ri_(a;)J+(-a;)) = 2^1im (e5^/"'CJ-(r,cu)) , (5.15d) 

(J+(a;)J_(-u;)) = 2 lim , (5.15e) 

(J_(a;)J+(-a;)) = 2 lim (e3^/^C-;-(r,cu)) , (5.15f) 

r^oo \ / 


with all other 2-point functions vanishing identically. In the x, y basis these become 

(T,A^)TU-^)) = {T„(w)n,{-w)) = 1 _,lim (C^lr, u,) + ""(r.i.,))) , (5.16a) 

(r,,(a.)T„(-a,)) = - (r,,(a.)T,,(-a<)) = (C+””(r,a,) - C^-"'‘(r,a.))) , (5,16b) 

(r,,(a.)J,(-a<)> = + C^rirM)) , (5,16c) 

(r„(a.)J,(-a.)) = - (r„(a.)J,(-a;)> = (e^c/l (C+;“(r,a,) - CJ”‘(r,a,))) , (5,16d) 

= (J,(w)J,[-^)) =i_,lim (c3-/5(cr“(r,a,)+C-"’‘('-,a.))), (5.16e) 

(M^)J,{w)) = - (J,{u,)U^)) =^_,lim (e*/5(cr“(r,a,)-C-"’‘('-,a.))). (5,161) 

The limits on the RHS can be evaluated explicitly by considering asymptotically AdS back¬ 
grounds of the form (3.1) so that 

a = ao + Eg i(Q - 2HUo)Le-^/^ + O , (5.17a) 

yl = I + O , (5.17b) 

/ = 1 - + O , (5.17c) 

and asymptotically S ~ Eg, Ft ~ Ho for some constants Eq and Flo. Note that this asymptotic 
form of / follows from the third equation in (3.6), even for backgrounds with running scalars. As 
we shall see, these asymptotic conditions for the background imply that the response functions 1Z± 
asymptotically behave as 


77±(r,a;) = 77±(i)(a;)e-^/^ + ©(e'^’’/^). 


(5.18) 
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where the functions TZ^(i'^{uj) are determined by the boundary conditions in the interior of the bulk 


spacetime. Evaluating the limits (5.16) using these asymptotics we hnally obtain 



{Ttx{uj)Ttx{-uj)) 


i—{Ttx{uj)Jyi-u})) -h {Ttt), 

UJ 

(5.19a) 


{Ttx{uj)Tty{-uj)) 

= ~ {Tty{io)Ttx{—Oj)) = 

i ^))) 

UJ 

(5.19b) 


{Ttx{ui)Jx{-i^)) 

= {'^tyi^)Jy{~<-^)) = 

i — {Jx{oj)Jy{-Lo)) - (j‘), 

UJ 

(5.19c) 


{Ttx{uj)Jy{-Uj)) 

= ~ {'^ty{^)Jx{—^)) = 

i {Jx{^)Jx{ ^))? 

UJ 

(5.19d) 

where 

Jxi ^)) 

= {Jy{ljj)Jy{—Uj)) = 


(5.20a) 


{Jxiuj)Jy{-U})) 


“ (^o7^Ci)(w) - 2nou;) , 

(5.20b) 

with 



T 

-H 

(5.21) 


and the renormalized 1-point functions {Tu) and (J*), corresponding respectively to the energy 
and charge densities, are given in (3.35) and (3.34). Notice that the only independent 2-point 
functions are the current-current ones, while all other non identically vanishing 2-point functions 
are expressed in terms of the current-current 2-point functions and the 1-point functions of the 
background. In fact, the relations (5.19) are nothing but the Fourier transform of the 2-point 
function Ward identities (2.41), and so they are purely kinematic. 

A number of comments are in order here. Firstly, note that the Ward identities (5.19) imply 
that certain 2-point functions can potentially diverge as cu —0 for non-zero magnetic field. However, 
as we shall see in Section 6.4, the leading behavior of the response functions 7^^(i)(u;) for small uj 
is (see (6.33)) 

^±(i)(^) = (2^no -Q)^ + (5.22) 

independently of the regularity condition imposed in the IR.® This ensures that the 2-point functions 
do not diverge in the small frequency limit for H ^ 0. Another remark is that all the expressions for 
the 2-point functions we have discussed so far, including the Ward identities (5.19), are independent 
of the IR boundary conditions. In order to compute the conductivities we will impose infalling 
boundary conditions on the horizon in order to obtain the retarded 2-point functions [71], but any 
other thermal 2-point function can be obtained by imposing the corresponding boundary condition 
on the horizon. As we shall show in Section 6.2, imposing infalling boundary conditions on the 
horizon implies that the response functions Tl±{io) are related via and hence 

7^+(l)(u;) =7^_(l)(-a;)^ (5.23) 

Finally, it should be stressed that we have computed the renormalized 2-point functions without 
ever evaluating the on-shell action, contrary to the usual procedure in the literature. What allowed 
us to do this is the fact that the renormalized canonical momenta are identified with the 1-point 
functions in the presence of sources via (2.37). This is practically very important in Lorentzian 
holography, since the proper evaluation of the on-shell action for thermal correlators involves mul¬ 
tiple boundaries [72-74] . 

®More precisely, it is independent of the IR boundary conditions, as long as the constant ci in Section 6.4 is 
non-zero. 
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5.2 Kubo formulas and transport coefficients 

The thermoelectric conductivities are defined in linear response theory via the relations 


( {J) \ ^ ^ T^\ ( E \ 

\{Q)I \Td.Tk)\-V\ogT ) ' 


(5.24) 


where J is the electromagnetic current and Q is the mixed thermoelectric current defined by 


{Qa) = {Tat) - a,b = x,y, (5.25) 

and fj, is the chemical potential. E here is the applied electric held and VT is the spatial temper¬ 
ature gradient. The linear responses in the currents are encoded in the matrix of thermoelectric 
conductivities, where 0 is electric conductivity, S and & are the thermoelectric coefficients and k 
is the thermal conductivity. Each of these quantities is a 2 x 2 antisymmetric matrix. The hats 
indicate that these quantities are the conductivities in the presence of a magnetization current [35] . 

Combining this dehnition of the conductivities with the expressions (5.13) for the 1-point 
functions of the stress tensor and electromagnetic current and the 2-point functions (5.16) we 
arrive at the Kubo formulas 

dab = (5.26a) 


Taab = iui ^{Qa{u})Jb{-‘^))R = ^{Tati‘^)Jb{-‘^))R- fJ^CFab, 


(5.26b) 


TKab = il^ ^{Qaiuj)Qb{-‘^))R = ^ {Qa{oj)Tu{-Oj))R - ytToiab 

= {Tat{oj)Tu{-uj))R - iJ,iu;~^{Ja{uj)Tbt{-uj))R - yTaab, 


(5.26c) 


where the subscript R indicates that the retarded correlators must be used. 

We can now use the Ward identities (5.19) -which we emphasize were derived holographically- 
to express the thermoelectric coefficients and the thermal conductivity in terms of the electric 
conductivities only. From (5.19) we get 


iujTOtxx — Tl^xy T P: 

UjToixy — 'i'R^xx ^P^xy-i 


(5.27a) 

(5.27b) 


where p = (J*) = —2QIk^ is the electric charge density we found in (3.32). Combining these two 
relations as 

OjT {cXxy T ^^xx) — xy T ^^xx) Ojp {^xy T ^^xx) T (5.28) 

and defining 

(7^ — ^xy T ’^^xxi — ^xy T ^^xx^ (5.29) 

we arrive at the expressions 

' ' (5.30) 


u)Ta± = —{ojp ± H)a± ± p, 


reproducing the result of [36]. However, here we have arrived at this result through a holographic 
derivation of the Ward identities (2.41), which hold in a generic theory including scalar operators. 
Similarly, combining (5.26c) and (5.19) we get 


ojTk± = —{ojp ± H)Ta± ± (e — pp), 


(5.31) 
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where R.± = Rxy =t iRxx-, and e is the energy density defined in (3.35). 

It follows that the only conductivities we actually need to compute are the electric conduc¬ 
tivities, which are holographically expressed in terms of the response functions through (5.20) as 

2i 

^ XX — ^yy — 2 (5.32a) 

Cfxy = -CTyx = ^ - 2110^ , (5.32b) 

where the response functions are computed with ingoing boundary conditions on the horizon. Note 
that there is no magnetization subtraction for the electric conductivities, which is why we have 
dropped the 

6 Response functions from the Riccati equation 

In the preceding sections we have shown that the renormalized 2-point functions and the cor¬ 
responding conductivities for generic asymptotically AdS backgrounds of the form (3.1) can be 
expressed in terms of the response functions 7^^(i)(w), which are as yet undetermined. These re¬ 
sponse functions are computed by solving the Riccati equations (4.15) and imposing appropriate 
regularity conditions in the interior of the bulk spacetime. Typically these equations can only be 
solved numerically, but certain analytic results can be obtained by taking various limits. In this 
section we determine the general solution of the Riccati equations in the small and large frequency 
limits, as well as the ultraviolet and infrared asymptotic solutions, before presenting an algorithm 
for obtaining the exact solution numerically. 

6.1 UV asymptotic solutions 

We begin by determining the UV behavior of the response functions TZ±(r,Lo). Inserting the UV 
expansions (5.17) for the background fields in the coefficients (4.10) of the Riccati equations we 
obtain 

gi{u:,H) = 1 + 0(6-'-/^), 5 o(^, ±R) = cu + ©(e'^^/^). (6.1) 

It then follows trivially from the Riccati equations (4.15) that the general UV behavior of the 
response functions takes the form 

7^±(r,a;) = + ©(e'^’’/^), (6.2) 

where 7^(1)^ (w) are the only integration constants of these first order equations. As we have seen, 
these integration constants, which are determined by imposing suitable regularity conditions in 
the IR, are the quantities encoding all the dynamical information in the 2-point functions and the 
corresponding transport coefficients. 

6.2 IR asymptotic solutions 

The IR behavior of the response functions TZ±{r,uj) depends crucially on the type of background 
considered. The ansatz (3.1) includes both zero temperature (no horizon) and finite temperature 
backgrounds and the results for the renormalized Green’s functions obtained in Section 5 are 
applicable to both types of backgrounds. However, the conductivities involve the retarded Green’s 
functions at finite temperature and so they are relevant observables for backgrounds with a horizon. 
We will therefore consider explicitly only the IR asymptotics for finite temperature backgrounds 
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of the form (3.1) here. Confining backgrounds of the form (3.1) at zero temperature are discussed 
e.g. in [22, 46]. 

Assuming the geometry exhibits a horizon at u = Uh (corresponding ot the smallest root of 
the equation f{u) = 0) giving rise to a finite temperature T, the function / increases from zero 
at u = Uh in the IR to the value / = 1 at u = 0 in the UV. In particular, f{u) admits a Taylor 
expansion near Uh of the form 

f{p) = AnTp + O {p^) , p = Uh-u, (6.3) 

while the warp factor satisfies = 0(1) as p —)• O"*". These asymptotic conditions, together 

with the assumption that S((/)) and n(y) remain hnite at the horizon, allows one to determine the 
leading asymptotic form of the coefficients (4.13), namely 

hi =--+ 0(1), ho = + O {p~^) . (6.4) 

It follows that the near horizon behavior of the general solution of the fluctuation equations (4.12) 

takes the form 

T±(p) = 4-(w)p-4^ (1 + 0(p)) + C?^\u:)p^ (1 + 0(p)), (6.5) 

where and are arbitrary integration constants, multiplying the two linearly independent 
solutions ~ respectively infalling and outgoing - of the second order equations (4.12). The retarded 
Green’s functions are computed by setting the outgoing mode to zero at the horizon, while the 
advanced Green’s functions correspond to setting the infalling mode to zero [71]. It follows that 
the near horizon behavior of the response functions 

7^± = drlog\£±\ = /^/4"^(9plog|T±|, (6.6) 


must be of the form 

n±{p,uj) 


—iuje "^(“'»)(47rrp) + C)(p^/^), Retarded, 

ia;e“"^(“'*)(47rTp)“^/^ + C>(p^/^), Advanced, 


(6.7) 


depending on whether we want to compute the retarded or advanced 2-point functions. These 
conditions on the horizon determine the sole integration constants 7^(i)^(a;) in the solution of the 
Riccati equations (4.15). Note that these boundary conditions are invariant under the combined 
transformation u —?■ —iv and complex conjugation. This leads to the relation (5.23) we mentioned 
earlier between the response functions. Another important remark is that the expansion (6.5) as 
written here is only strictly valid for nonzero ui. This is because the 0{p) terms that appear in this 
expansion contain inverse powers of ui, thus rendering the limit w —)• 0 ill defined. We need to keep 
this fact in mind when determining the small cu behavior in Section (6.4). 


6.3 Universal large ui solution 

Besides the UV and IR asymptotic expansions, the large and small frequency solutions of the 
fluctuation equations can be obtained analytically. These determine respectively the large and 
small frequency behavior of the response functions 77.(i)^(a;), and hence of the conductivities. 

In the large frequency limit the coefficients (4.13) become 

ho = ‘^ + 0{io), 4 = 41 og|S/| + 0(a;-2). (6.8) 
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Assuming these expansions hold uniformly in [0,Uh\, the fluctuation equations (4.12) to leading 
order in uj simplify to 

=0, (6.9) 

whose general solution is 

£± = {ci + 0{uj~^)) exp (^iuj + {c 2 + 0{u}~^)) exp tImo) ’ 

with Cl and C 2 arbitrary integration constants. From the near horizon behavior of / in (6.3) we 
deduce that infalling boundary conditions on the horizon corresponds to setting C 2 = 0. This yields 


= 9,.log|T±| =^9„log|T±| =-icj/ + C>(a;°), 


and so from (6.2) we conclude that 


’^±(i)(‘^) — 


—iu! + 0(uj^), Retarded, 
iuj + 0{uj^), Advanced. 


( 6 . 11 ) 


( 6 . 12 ) 


From (5.32) then follows that for large frequencies the conductivities behave as 




" 1 

2^^ 0(cC ), (Txy^OJ^ — O'yx{oj') 



Uo + Oiuj-^). 


(6.13) 


Notice that this result, which agrees with e.g. eq. (6) in [36] and eq. (67) in [30] in the respective 
conventions, is independent of the particular background considered here and is thus universally 
valid for asymptotically locally AdS backgrounds of the form (3.1), generalizing previous results by 
adding potentially running scalars with arbitrary potentials. Moreover, as expected, it is indepen¬ 
dent of the magnetic field and charge density, which are relevant deformations of the theory and 
hence do not affect the UV physics. 


6.4 Universal small lo solution and its Fade approximant 


A universal result for the response functions 7^^(i)(a;) can also be obtained in the small frequency 
limit. From (4.13a) and (4.13b), using the background equations (3.6), we deduce that 


1 /p2A\' 

hi{uj, H) = 2A' -b V/S j ^ ’ 


(6.14a) 


^-2A 


ho{uj, ±H) = T ^ 


-2A /„2A„/\' 


e a 


f 


LO 


w-b ^=F 


la'/ 


f2 ' 4^3 f 


(2A \ ' 

^—jco^ + 0{A). (6.14b) 


Inserting an expansion of the form 

= 4°^+ w 4^^+o {A) , 


(6.15) 


- 30 - 










in the fluctuation equations (4.12) leads to the following equations up to 


0{u^) 

0{u^) 

0(a;2) 

0(a;3) 


—24 —24 

^.( 0 ) ^ _ e_ ^ 


2A f2 I 


( 0 )' 


' 




= 0 , 


( 1 )' 




l\ ' 


H 


2A„/\ ' 


f fe^^a‘ 

= ±4' 


/ 


s. 


(0) 
± > 


( 2 )' 


^2Af2 ( ^ 


' 


.^(0) 


o2A , 


= —- 


f / J2A\> , , .f / J2A^i\' 

J I ^ \ c' 0) I J / e a \ 


4^2 V/s y 


\ c'(0) , 

I 4 


/ 




( 1 ) 
± ) 


(3)- 


^2Af2 ( ^ 


/\ ' 


^2A 


= ± 


rO! 


4/^3 V s/ y 




2A . 


„2A\' , , „2A , , o2v4 f /„2A\' , , . 

I ff" - f?> ± 4 


/ 


4 i 42 V /s y 


(6.16a) 

(6.16b) 

(6.16c) 


/ ye^^a'V (2) 

■ 


H 


f 

(6.16d) 


Note that the 0{u!^) equation is a homogeneous second order equation and hence the corresponding 
general solution contains two integrations constants. Moreover, the subleading in uj equations are 
second order inhomogeneous equations but with the same homogeneous solutions as the 0{uj^) 
equations. Without loss of generality we will absorb all freedom in choosing the homogeneous 
solutions into the 0{uj^) solution by allowing the integration constants to potentially depend on 
the frequency. 

At 0{uj^) the two linearly independent solutions are 


^ =Cif + C2fp 


, p{u) = r 

Jo 


dn' 


(6.17) 


where the lower limit of integration in p{u) has been chosen so that the integral is well defined. In 
the UV these behave as 


f{u) = 1 - m4 + O {u'^) , p{u) = — + 0 (4) , 
while near the horizon we have 


— 2A{uh) 

f{u) = AttTp + O {p'^) , p{p) = +0{\ogp). 


(6.18) 


(6.19) 


Both / and fp are therefore regular on the horizon and so we cannot a priori exclude any of the 
two solutions at this order in the frequency. As we shall see, we need to determine the expansion 
(6.15) up to O{oj‘^) in order to find suitable linear combination of the integration constants ci and 
C 2 that corresponds to the desired IR boundary conditions. The inhomogeneous solution at O(w^) 
takes the form 


c(i)_l_ 


du 

e2^/2 


Cl 
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dM + C2 / fp 


f 




f 


du 


( 6 . 20 ) 


where again the lower limits of integration have been chosen so that the integrals are well defined. 
Near the horizon these the 0{ui^) solution behaves as 


yi) 

•'ll 


(Q - 2HIlh) / 


HT. 




^-2A{uh) \ 

Cipl0gp-C2^^^^l0gpl , 


( 6 . 21 ) 
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where II/j = n(x(ri/i)). In the UV we have 

4'^ = Cl {Q - 2^no) u + O + C20(u^). (6.22) 

Combining the O{io^) and O(uj^) results, near the horizon the small frequency expansion behaves 

^-2A(uh) 


as 


S± ~ Cl [AirTp + ujO{plogp) + 0 (w^)) + C 2 


AttT 


+ ujO{logp) + 0{uj 


(6.23) 


It follows that for any value of the integration constants ci and C 2 this expansion breaks down when 
— |a;| log /9 « 1. However, as long as p » the expansion (6.15) is well dehned. Going one 

order higher, the 0{u‘^) inhomogeneous solution takes the form 


4'^ = -/ f 

Jo 


du 


( 0 ) , ,( 0 ) / 


f 


" “ ' 4 ) ) dh, 


(6.24) 


where the lower limit of the first integration, 0 < u* < u/i, is an arbitrary point between the horizon 
and the boundary since the integrand diverges both at the horizon and the boundary. Near the 
horizon this behaves as 


4^ = ci(-^ 


^2A{uh) 


where S/j = S((()(n/i)), while in the UV 

( 2 ) ^ 3L^ 

± 4F2So 


log p +0(1) + C20((logp)" 


Mcl - ) U + 0{u^). 


(6.25) 


(6.26) 


At this point we can determine the relation between the integration constants ci and C 2 cor¬ 
responding to ingoing boundary conditions on the horizon. Provided p » the ingoing 

solution in (6.5) can be expanded as 


£± =d±{uj) —u!logp + 0{u!logpf) 

On the other hand, in the same limit the small u solution we found behaves as 

g-2A{ufi) 


S± ~ C2- 


dvrT 

and hence we conclude that 


(Q-2HUh), Cl 


Cl = C2- 


(47rT)2t9 


iHJ^h V (Q - 2HUh) . 


(6.27) 


(6.28) 


(6.29) 


( 2 ) 

Applying the same argument to subleading orders in the expansion determines the point n* in <51^ k 
Finally, observing that 

4'^ T 1^4^ = Oiu% (6.30) 

as It —>■ 0, it is straightforward to show that 4^^ = O(u^) and hence in the UV we get 

1 C 2 ' 


= Cl 1 ± 


CO 


H'Lf 


[q 


— IHIin U + 


“iL^co^ 

4772Sr 


M - 


3L2 


Cl 


u + 0{u'^,cj‘^) 


(6.31) 
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It follows that 


7^± = 9.log |^±| = -f/^e-^dulog\£±\ ~ I (g - 2HUo) 

from which we read off the response functions 


3-L^u;^ 

4F2So 




3L2 Cl y y ’ 
(6.32) 




3ML‘^u‘^ (47rr)2e4^(“^)w3 ^ 

4/72So + 16So//3(iFSi, ^ (g - 2^^) ^ 


(6.33) 


Notice that the leading O{oj) part of this expression is universal and independent of the IR bound¬ 
ary conditions as long as ci 7 ^ 0. This temperature independent part of the response functions 
determines the universal DC Hall conductivity and it is crucial to remove the apparent pole at 
cj = 0 in the conductivities (5.32). From (5.32) it follows that the small frequency behavior of the 
conductivities is 


_ 3ML2 

<yXX - Oyy - 

_ P__ 

O'xy — O'yx — 


(47rr)2e4^(“'^)S;,u;2 

iuj —--T (D ioj ) , 

+ {Q- 2HUh)^) 

(47rr)2e4^K)(g-2Hnga;2 ^ ^ 

8 H3 k2 (H'2S2 + (g - 2Hn;,)2) ^ ’ ■ 


(6.34a) 

(6.34b) 


Given that M is related to the temperature T, the second order correction in w of the response 
function brings in a temperature dependence in the leading nontrivial correction of the transport 
coefficients. From the Ward identities (5.19) then we obtain 


{Ttx{uj)Jx{-(^)) = {Tty{uj)Jy{-Uj)) = 


{Ttx{uj) Jy {-Oj)) = -{Tty{uj)Jx{-Uj)) = 


( 4 ^r) 2 ^ 4 AK)(g_ 2 ffnga;^ ^ ^ .3^ 

+ {Q- 2HUhy) ^ ^ ’ 

3ML2 (47rr)2e^^(“'*)S/,a;2 

2Hk^ ^ 8Hk2(h2s 2 + (g - 2H'ng2) 


(6.35a) 

+ O , 

(6.35b) 


{Ttx{uj)Ttx{-uj)) = {Ttyiio)Tty{-io)) = -e + -^ + 


(47rr)2e^^(“'‘)S,jia; 


{Ttx{uj)Tty{-u;)) = -{Tty{aj)Ttx{-uj)) = 


2^2 8k2(h'2s 2 + (g - 2Hng2) 

(47rT)^e^-^(^'-)(g - 2Hnfe)^a; 

8Hk2(h2s 2 + (g - 2Hn,,)2) 


+ (cu^) , 

(6.35c) 

+ O (u;2) . (6.35d) 


These agree with equations (49)-(51) in [34] for the dyonic Reissner-Nordstrom black hole, but the 
present derivation holds even for backgrounds with non-trivial scalar profiles. 


Fade approximant 

The small frequency expansion (6.33) of the response functions can be considerably improved by 
means of a Fade approximant, which is capable of capturing certain poles of the conductivities. 
It turns out there are two different Fade approximants that correctly capture the behavior of the 
conductivities at different limits of parameter space. The two Fade approximants correspond to 
two different terms dominating the response function (6.33). Writing (6.33) as 

7^(l)±(a;) = + 0(a;^), (6.36) 
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the two relevant Fade approximants are 








.[ 1 ] 


7^ 


[ 2 ] 


a ;2 » 


'^\i)±^ + 


,, 2 ('- p [21 >,2 

,-7^® 




(6.37) 


''{i)± '''(i)±“ 

The first Fade approximant is a good approximation to the response functions at low temperature, 
i.e. near extremality, and leads to a pole at 


^*dl 


4g(Q - 2gno) 
3L2M 


(6.38) 


which agrees with (81) of [36] in the small H limit keeping w/i^ fixed. The second Fade approximant 
leads to poles at 




12 LUfg 

(47rT)2 


{iHJ:hT{Q-2HUh)), 


(6.39) 


which agrees with (74) of [36] in the hydrodynamic limit where H and Q are sent to zero keeping 
H/uj and Q/oj fixed. As we shall see in Section 7 these Fade approximants capture the poles closest 
to the origin of the complex oj plane to very good approximation. 


6.5 Numerical solution of the Riccati equation 

The Riccati equations (4.15), or equivalently the second order linear equations (4.12), are not in- 
tegrable in general and so one needs to solve these equations numerically. One can solve either the 
Riccati equations or the second order equations, but the fact that the Riccati equations directly 
determine the response function by imposing only IR boundary conditions is a clear advantage 
compared to the second order equations where one must keep track of the arbitrary source.® More¬ 
over, the relation 77+(w) = 77_(—w)* implies that instead of computing 77+ and 77_ for positive uj 
we can compute only 77+(a;) for all w. Thus, in order to compute the conductivities one only needs 
to solve one first order ordinary differential equation which is simpler than solving a set of coupled 
second order differential equations as was previously done in the literature. 

In order to solve the Riccati equations (4.15) numerically it is convenient to introduce the new 
dependent variables 

(6.40) 


0+ = 0-1 ±uj{Q- 2IiH^ 


in terms of which the Riccati equations take the simpler form 


/ (770+ + - ^77^ - S”! 

W0+ T 

'q - 2UH^ 

)A„, 


For retarded Green’s functions the horizon condition (6.7) translates to 


(6.41) 


0+ ~ w-i ±{Q- 2IihH) ) , (6.42) 

on the horizon, while in the UV 

0+ ~ w-2 (^RSo77(i)+ ±u{Q- 2noR) ) , (6.43) 

®It is worth pointing out that the boundary value problem for the second order equations (4.12) is well defined, 
despite the apparent singularity of the the coefficients ho and h\ at the zeros of F = . Indeed, it is 

straightforward to see that ho and hi have a simple pole at the zeros of F and therefore these correspond to regular 
singular points of the second order equations [22, 75]. 
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from which one can read off the renormalized response functions Equation (6.41) can be 

integrated using any standard solver for ordinary differential equations. However, the horizon is 
not a regular singular point and imposing boundary conditions on the horizon requires some care 
in the numerical analysis. A standard technique is to use a Taylor expansion in the vicinity of the 
horizon, and match the numerical solution at some small distance away from the horizon. However, 
we found that even a near-horizon expansion to 0{uh — u)^ was not sufficient for stabilizing the 
numerics for (6.41) with NDSolve in Mathematica.^® Instead, using a Fade approximant based on 
the near-horizon expansion to 0{uh — u)^ in order to impose the IR boundary condition sufficiently 
away from the horizon worked very well with NDSolve. 


7 Example: Dyonic Reissner-Nordstrom black hole 


In this section we will apply the above general analysis to the dyonic Reissner-Nordstrom black hole 
(3.17), which was first studied in [34, 36]. The dyonic Reissner-Nordstrom black hole is a solution of 
the background equations (3.6) either for a model with constant potentials, i.e. V{4>,x) = —6/L^, 
S(</)) = Eo, ZicP) = Zo and n(x) = Hq as in [34, 36], or in a theory with arbitrary potentials 
that admits asymptotically AdS solutions, provided the scalars are set to their vacuum AdS value, 
which can be taken without loss of generality to be zero. Setting the scalars to zero is a consistent 
truncation of the background equations (3.6) provided V{4>,x), Z,{4>) and n(x) don’t have a linear 
term in their respective variables in a Taylor expansion around 4> = x = 0. However, while in the 
theory with constant potentials the dyonic Reissner-Nordstrom black hole - as we shall demonstrate 
momentarily - is the only asymptotically AdS solution of the form (3.1), this is not generically the 
case in the theory with non constant potentials. In particular, a generic theory will admit in 
addition to the dyonic Reissner-Nordstrom black hole a hairy black hole with the same charges and 
a phase transition between the two solutions will generically occur at some critical temperature. 

In order to show that the dyonic Reissner-Nordstrom black hole is the only asymptotically 
AdS solution of the form (3.1) in the case of constant potentials, let us consider the background 
equations (3.6), which in this case reduce to 


(e2^/H')'- ^ + So (R2 + a'2) = 0, (7.1a) 

2 (H" - A'2)+</,'! +Zox'l = 0, (7.1b) 

(e^^/)' - 4So (R^ + a^) = 0, (7.1c) 

= 0, (7.1d) 

(e'^/Xs)' = 0, (7.1e) 

(Soa'+ 2noR)'= 0. (7.1f) 

The general solution of the last equation is 

a{u) = ao — Qu, (7.2) 


where the constant Q is related to Q that we introduced before by SqQ = 2noi7 — Q- Moreover, 
the scalar fields can be expressed in terms of f{u) and A{u) as, 


(f>B{u) = (po + C<l> / 
Jo 


du' 


g2A{u') 


Xb{u) = Xo + C^ 

Jo 


du' 




(7.3) 


also solved (6.41) using the ODE integrator from the Python library scipy, finding the same results. 
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where (/>o, XO) C(j> and are integration constants, while using the solution for a, the blackening 
factor can be expressed in terms of A{u) as 

f{u) = h+ r du' (3M + 4So + Q^) n') (7.4) 

40 

where M is another integration constant. Finally, eliminating / from the first two equations we 
obtain a second order equation for the warp factor A(u), which is the only remaining unknown 
quantity, namely 

{Cl + ZoCl) {A" + 2A'^)^ + 2 {A" - A'^) (^3MA' + {4A'u + l)^o{H^ + Q^) - = 0. (7.5) 

However, since we are interested in asymptotically AdS solutions we must have A ~ log(L/ri) as 
u —)■ 0. The last equation then requires that = 0 for asymptotically AdS solutions and 

hence the scalars must be necessarily constant and so, without loss of generality, we can set them to 
zero. The resulting solution is the dyonic Reissner-Nordstrom black hole (3.17) in four dimensions 
for which 

(Pb{u) = 0, A{u) = log{L/u), f{u) = l-Mu^ + ^{H^ + Q‘^)u‘^, 

Xb{u) = 0, a{u) = ao - Qu. (7.6) 


The horizon radius is the smallest positive root of the quartic equation f{uh) = 0, i.e. 

Mul = l + ^{H^ + Q^)ul 


(7.7) 


This expression can be used to express the mass M in terms of the horizon radius Uh, namely 


<“> = - 5) - § (i - 


(7.8) 


Finally, the Hawking temperature can be found by the usual argument demanding that the Eu¬ 
clidean section be free of conical defects, leading to the expression 




(7.9) 


while the energy and pressure densities, respectively (3.35) and (3.36), become 

(7. 

In Figures 1 and 6 we plot the real and imaginary parts of the response function 77_(i)(t(j) as a 
function of a real frequency for two different choices of the electric and magnetic fields, correspond¬ 
ing to different temperatures. The values of the electric and magnetic fields in Fig. 1 correspond to 
a nearly extremal black hole, where the first Fade approximant in (6.37) provides a good approx¬ 
imation of the response functions for small frequencies. In particular, the Fade approximant very 
accurately reproduces the pole closest to the origin of the complex frequency plane as can be seen 
explicitly in the plot. The conductivities for these values of the electric and magnetic fields are 
plotted in Figures 2, 3 and 4. The fact that there is no Drude peak and the poles are not located 
at zero frequency is a consequence of the broken translational symmetry due to the magnetic field 




e = 


V = — 


ML‘^ 


e + 2V = Q. 
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Figure 1. Plots of the real and imaginary parts of the response function as a function of a real 

frequency w for i/ = 0.6, Q = 1.6156, T = 0.00239, together with the small frequency expansion (6.33), 
the first Fade approximant in (6.37), as well as the asymptotic behavior (6.12). Both the location of the 
horizon Uh and the AdS radius L are set to 1. As is evident from this plot, the Fade approximant is a drastic 
improvement of the small w expansion in this regime of parameter space. In particular, it captures very well 
the pole in the Real part of 7?._(i)(a;), but not the delta function at the same frequency (as predicted by the 
Kramers-Kronig relations) in the imaginary part. 


[34, 36]. How the Drude peak is recovered in the limit of vanishing magnetic field is illustrated in 
Fig. 5. The values of the electric and magnetic fields in Fig. 6 are instead in the hydrodynamic 
regime where the second Fade approximant in (6.37) provides a good approximation of the response 
functions for small frequencies, including the location of the poles, which have now moved away 
from the real axis according to (6.39). The corresponding conductivities are plotted in Figures 7, 
8 and 9. 

In Fig. 10 we compare the location of the pole of the conductivities Gxx and CTxy nearest 
to zero on the real frequency axis as a function of the magnetic field FT at a number of different 
temperatures with the result (6.38) predicted by the first Fade approximant in (6.37). As expected, 
the Fade correctly gives the location of the pole at very small magnetic field, independently of the 
temperature. However, the lower the temperature the agreement extends to a higher value of the 
magnetic field. An analogous plot for the poles (6.39) in the hydrodynamic regime can be found in 
Fig. 3 of [36]. 

Finally, in Fig. 11 we plot ju+j as a function of the complexified frequency for a number of 
values of the electric and magnetic fields keeping + Q‘^ = 1 fixed, reproducing Fig. 1 in [36] 
(see also Fig. 9 in [26]). To capture the non-liner regime where the cyclotron poles deviate from 
the semi-circle configuration it was necessary to impose the IR boundary condition using a Fade 
approximant in the near horizon expansion, as discussed in Section 6.5. This transition regime is 
not visible in the plots of [36] and [26] . 
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Figure 2. Plots of the real and imaginary parts of the conductivities and Uxy as a function of a real 
frequency oj iov H — 0.6, Q — 1.6156, T — 0.00239, together with the small frequency expansion (6.33), the 
first Fade approximant in (6.37), as well as the large oj asymptotic behavior (6.13). Again the location of 
the horizon uu and the AdS radius L are set to 1. The location of the poles in Imaxx and Reaxy are well 
approximated by the Fade approximant and are given by (6.38). However, the delta functions in Recra;^; and 
Imaxy, which are related to the aforementioned poles via the Kramers-Kroning relations, are not captured 
by the Fade approximant. As —>■ 0 the delta functions move towards the origin of the complex to plane, 
giving rise to the well known Drude peak, which reflects translation invariance. The fact that the poles (and 
hence the delta functions) in these plots are away from the origin is a consequence of broken translation 
invariance due to the magnetic field [34, 36]. A similar effect occurs when translation invariance is broken 
by impurities [76]. 

8 Concluding remarks 

In this paper we presented a general framework for the holographic analysis of asymptotically AdS 
backgrounds with finite charge density and a constant magnetic field, including the systematic com¬ 
putation of the renormalized 1- and 2-point functions and the corresponding transport coefficients. 
The importance of such dyonic backgrounds in the holographic study of strongly coupled systems 
in both condensed matter and high energy physics was our main motivation for carrying out the 
analysis in a general and systematic fashion, in the hope that our results can be directly used in 
numerous applications. 

There are three important aspects in our general holographic prescription that we have tried 
to emphasize, all of which rely on a radial Hamiltonian formulation of the bulk dynamics. The first 
is a general recursive prescription for holographically renormalizing the theory. This is particularly 
important in the presence of running scalars, which can contribute to the UV divergences. More¬ 
over, correctly renormalizing the theory is crucial to ensure that ultralocal and quasilocal terms in 
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Figure 3. Plot of |cra;a;P as a function of a complex a; for = 0.6, Q = 1.6156, T — 0.00239 (left), together 
with the first Fade approximant in (6.37) (right). As is evident from these plots, the Fade approximant 
captures the poles closest to the origin extremely accurately. 



Figure 4. Plot of \<7xy\‘^ as a function of a complex lu foi H = 0.6, Q = 1.6156, T = 0.00239 (left), together 
with the first Fade approximant in (6.37) (right). 


correlation functions (and hence transport coefficients) are compatible with the Ward identities. 

The second aspect we wanted to highlight is the fact that correlation functions are much more 
efficiently extracted directly from the solutions of the equations of motion instead of first evaluating 
the on-shell action and then taking derivatives. The (renormalized) radial canonical momenta are 
holographically identified with the 1-point functions of the dual operators in the presence of sources, 
which can be further differentiated with respect to the sources to give any desired n-point function. 
The only instance one actually needs to compute the on-shell action is in the computation of the 
free energy of the background solutions. 

The third aspect of our analysis that we consider important concerns the fluctuation equations 
for determining the 2-point (and higher) functions. These are in general a system of coupled 
linear second order equations and require boundary conditions in the UV and regularity conditions 
in the IR. Unless one can decouple these equations it is a priori tricky to identify which modes 
are the independent sources and which the responses, although in principle this can be addressed 
systematically using the symplectic form of the bulk theory [44] . However, there is a straightforward 
alternative which utilizes the manifest symplectic structure of the Hamiltonian formalism. Namely, 


- 39 - 

















UJ UJ 


Figure 5. The parameters used in these plots are Q = Ij -^ = 3, So = 1; IIq = 0 and L = 1. Left panel: 
Tieaxx as a function of w for several values of the magnetic field H. As H decreases, the real part of axx 
tends to behave as a delta function keeping the area below the curve fixed but rising sharply and becoming 
steeper. This behavior occurs for small but nonzero w but for even smaller lo it vanishes steeply. Right 
panel: Im cra;^; as a function of oj for several values of the magnetic field H. As H decreases, the imaginary 
part of (Jxx tends to behave as ^ S for small but nonzero uj but for even smaller uj it vanishes steeply. We 
conclude that the magnetic field regulates the delta function which is present at w = 0 for zero magnetic 
field in which case there is a Drude peak involved. In particular, the Drude peak predicts Imaxx 1/w and 
RecTajj, ^ <5(a;) which is what we find by taking a sequence of decreasing H configurations. 




Figure 6. Plots of the real and imaginary parts of the response function 7?._(i)(w) as a function of a real 
frequency w for iJ = 0.28, Q = 0.04, T = 0.2324 and Uh = L = 1, together with the small frequency expansion 
(6.33), the second Fade approximant in (6.37), as well as the asymptotic behavior (6.12). Again, the Fade 
approximant is a drastic improvement compared to the small w expansion in this regime of parameter space. 
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Figure 7. Plots of the real and imaginary parts of the conductivities and a^y as a function of a real 
frequency w for = 0.28, Q = 0.04, T = 0.2324 and Uh = L = 1, together with the small frequency 
expansion (6.33), the second Fade approximant in (6.37), as well as the large uj asymptotic behavior (6.13). 
Again, the location of the peaks is very well approximated by the Fade approximant and is given by (6.39). 
Note that the poles have moved away from the real axis, which is why only peaks appear in these plots. 


by trading the linear second order fluctuation equations for the corresponding first order Riccati 
equations one automatically eliminates the sources out of the problem and computes directly the 
correct response functions by imposing only regularity conditions in the IR. Besides automatically 
taking care of the identification of the sources and response functions, the Riccati equations can be 
used to directly holographically renormalize the 2-point functions (only computing the terms that 
are contributing to the particular 2-point function [42]), and they render the numerical solution of 
the fluctuation equations considerably simpler by eliminating the arbitrary sources from the start 
[43]. 

Finally, the radial Hamiltonian formulation of the bulk dynamics often leads to powerful so¬ 
lution generating techniques for background solutions, as we have demonstrated in Section 3.2 
and Appendix A. In particular, we found infinite families of exact RG flows interpolating between 
AdS in the UV and hyperscaling violating geometries in the IR, some of which we expect exhibit a 
gapped and discrete spectrum of fluctuations. Interestingly, even the purely electric version of these 
solutions is not strictly included in the classification of hyperscaling violating geometries discussed 
in [61, 68, 70] since those assume a single exponential behavior for the scalar potential. However, 
the exact solutions we present here involve subleading terms in the scalar potential in an essential 
way. 

Computing the conductivities in these backgrounds is one of our immediate priorities [46]. 
Another interesting question is whether some of these solutions can be embedded in gauge super- 
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with the second Fade approximant in (6.37) (right). The plots on the top are zooming in on the poles closest 
to zero visible in the plots at the bottom. Again, the Fade approximant captures the poles closest to the 
origin extremely accurately. 

gravity. The potential V{(j)) can be easily embedded in supergravity by choosing Wo{(l)) to be (for 
example) the supersymmetric superpotential. However, the gauge kinetic function 'S(<p) is then 
determined as well, which renders the embedding not a completely trivial question. Finding exact 
families of dyonic black hole solutions with scalar hair and adding momentum relaxation are two 
other interesting directions we plan to explore [77]. 
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Figure 9. Plot of lo'xyl'^ as a function of a complex lu for H = 0.28, Q = 0.04, T = 0.2324 (left), together 
with the second Fade approximant in (6.37) (right). The plots on the top are zooming in on the poles closest 
to zero visible in the plots at the bottom. Again, the Fade approximant captures the poles closest to the 
origin extremely accurately. 

Appendices 

A Exact superpotentials and RG flows in various dimensions 

In this appendix we provide a few more examples of solutions to the superpotential equation (3.14), 
leading to exact RG flows in various dimensions. 

(i) d = 3: 

The solution (3.18) can be generalized to include a non-trivial axion, namely 

W{A,cP,x) = Wo{(P,x)V^T7^, (A.l) 

with the potentials given by 

V{cj>,x) = + Z-\cl>)Wox - 2^o^ 

R2s(</,) + (g - 2n(x)R)2E-i(</>) = ^ (w^o^ + Z-^mo^ + • (A-2) 

However, contrary to the purely dilatonic case discussed in Section 3.2, the function lTo(</>, x) 
cannot be chosen arbitrarily here since it is constrained by the second equation in (A.2). 
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Figure 10. Plot of the location of the pole of Uxx and dxy closest to zero on the positive real w axis as a 
function of H (solid lines), compared with the pole (6.38) predicted by the first Fade approximant in (6.37). 
As expected, the agreement is best at small magnetic field and small temperatures. 
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Figure 11. Density plots of |ct+| as a function of complex frequency for = 0 and Q = 1 (left), H = Q = 
l/v^ (center), and H = 1 and <5 = 0 (right) at T = 1I2 tt (cf. Fig. 1 in [36] and Fig. 9 in [26]). The 
white regions correspond to poles, while the blue regions to zeros. The location of the poles closest to zero 
forms a semicircle, along which the poles move as the values of H and Q are shifted, keeping T fixed [36]. A 
3-dimensional version of the plot in the center is shown in Fig. 9. 


(ii) d arbitrary, H = 0, Q 0: 

For zero magnetic field the superpotential equation (3.14) admits a solution of the form 




in any dimension, with ^ = q^L ^ and 

y(0) = IFf - E(,/.) = 


2Sn 




(A.3) 


(A.4) 
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Figure 12. Plots of jcr^xl^ (left) and \(Jxy\^ as functions of complexified frequency ior H = Q = 11^/2, 
T = 1/2 t:. Note that the poles are arranged in a configuration between the near extremal case in Figures 3 
and 4 and the hydrodynamic limit in Figures 8 and 9. 


As in the example in Section 3.2 the function Wo{4>) can be prescribed at will. The corre¬ 
sponding metrics again flow to a hyperscaling violating geometry in the IR with exponents 

2d — 3 

l<z<——6 = d + z-2, (A.5) 

and it is important to keep in mind that d here includes time. 

(iii) d = 5, Q — 2IioH ^ 0: 

In this case a possible solution of (3.14) is 

iy(A, <t>) = Woicj)) + {Q- 2UoH)Wi{(^)e-^^, (A.6) 


where 


VicP) = Wlf --Wo^ 




ITf + 


+ 5 ^ 2 ^ (^W^Wl + ^WoWi'^ = 


2H^ 


Q - 2noH 


(A.7) 


The last equation should be viewed as a differential equation for and the ratio j (Q — 

2noi7) must be treated as a parameter of the theory, unless it drops out of the combination 
. Again, Wo(</’) can be any function. 


B Gauss-Codazzi equations 

In this appendix we write for completeness the full set of Gauss-Codazzi equations for the model 
described by the action (1.1) in the gauge 

ds^ = dr^ -h dx^dx^^ A = Aidx'', (B-1) 

which is used throughout our analysis. These equations are the starting point for deriving the 
general fluctuation equations, which we give in Appendix C. The Gauss-Codazzi equations following 
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from the equations of motion (1.4) are 

Einstein rr: 

= i?[7] + 02 _ ^ z{cP) (x^ - d\dix) 

+ S(0) - E(0, x), (B.2a) 

Einstein ri: 

D,Kl - D,K = + Z{c^)xdiX + 2S(0)F,,7^''=ifc, (B.2b) 

Einstein ij: 

k] + KK] = i?;.[7] - - Z{^)kxdjX - 2E(0) (j^^AkAj + 

+ - E(0, x)) , (B.2c) 

Scalar 0: 

0 + K0 + 00 - ^ (x^ + 0'x5ix) + ( 27 '^'iiii + ) = 0, (B.2d) 

Scalar x- 

Z{x + Kx + Dx) + {x^ + 5V0.x) - ^ (^x + ^'nxF^'^kFjk) = 0, (B.2e) 

Maxwell i: 

^dr [V=X (S(0)7'^'ii + n(x)e*^'=E,fc)) = D, (s(0)F*^- + 2n(x)6*^’"ifc) , (B.2f) 

Maxwell r: 

A (S(0)7''ii + n(x)e'^''^Afc) = 0. (B.2g) 

C General fluctuation equations for d = 3 


Although in the main body of the paper we consider only certain time-dependent fluctuations 
around the backgrounds (3.1) with no spatial dependence, in this appendix we provide the complete 
set of fluctuation equations for d = 3 following from the Gauss-Codazzi equations in Appendix B, 
with generic fluctuations around the backgrounds (3.1) of the form 

lij ='yFj{k + hij{r,:>A, Ai = Af (r,x) + ai(r,x), 0 = 0B(r) + 9 ?(r, x), x = ^^(O + "^(a x). 

(C.l) 
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Introducing the quantities Sj = S = SI and S_\_ = SI — Sq, we have the following useful 

identities; 

-Rjjbs] = 0, Df = di, (C.2a) 

■ / . (0) f (0) 1 (0) . / 1 / \ \ 

K!=A5i + ^6l6^, K = dA + ^, K,, = -^,,=e‘^^(^A6,j--[f + 2A{l + f))5,oSjoj, 

(C.2b) 

( 0 ) . . 

(7fci7"') = (4i' + 2i)<5^' + (C.2c) 

1 . f . . (1) 1 (1) 

Ki = ^Si + ^{Sf6l-Sl6^), K = -S, R=d^djSi-aBS, (C.2d) 

= 2 {djSi + dkSi - d%k) , Rl = ^ (d^diSi + d^d^S^ - ObSI - d^d^S) . (C.2e) 

Expanding the Gauss-Codazzi equations in Appendix B and using these identities leads to the 
following set of fluctuation equations for d = 3: 

Einstein ij: 

(52 + l^dA + 1 /-v) dr + Db) si - f-^f [s^sl - SiSf) + ^S6^5i - [d^dkSf + d,d^si - d^d,s) 
+ 4Se-^^ {e^^r^oA + H^) {s^5i - 5^5°) + dj (^AS + + V^r)^ 

+ 4S^(^ + 4S - dd°5^) + -fi^d^adk) 

- 4idEe-2^ (Fe-2^ ((d^df + d^df)(5," + Sy) + di{S^Jf + ^^df)) 

-l^B (dpai - dia,)) 

- ^d/Se-2^ {-2f-^dta + r^S^a^ + 2He-^^{d^ay - dya^) - e-^^H\S^ + S^)) 

- jXe-2A^iE^<f {-r^d^ + H^e-^^) = 0, (C.3a) 

Einstein ri: 

dj (^i + f-^fisfsi - 5^d0) + Ir^fd^d^s'^ -d,(^s+ Ir^fS^o + 2 <Pb^ + 2Zxbt^ = 

- 4Ee-^^f-^d(diat - doaA + 4J:He-^^{5fdy - dfd^) - 4Si7e-2^(df 5° - df5°)d, (C.3b) 

Einstein rr: 

(d - i)As + ^r^fS± = {d^dj - dpB) si + 2P^ + 2ZxBf + {zpl -V^)ip- V^T 

- 2p^e-^^ (/"^d^ + H^e-p 

+ 2e-2^E {r^d^S^ - 2/-^ddt - 2He-^^{d^ay - dya,) + H^e-^^{S^ + Sp , (C.3c) 
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Scalar <y9: 

(p + 2{dA + - Z^XBT + ^^bS - ^ + Z^^xb) ^ - \^x<P'^ 

- {H^e-^^ - r^a^) 

- S^e-2^ {-2f-^adt + - H^e-^^{S% + S^) + 2He-^^{d,ay - dya,)) = 0, (C.3d) 

Scalar r: 

Z ^ + {dA + -/ ^f + Z + Obt^ + Z^xb'P + -^ZxbS — 2^xx'^ 

+ (^z^ 

+ 2 /-^/ 2 g-dA (_2n^^rdif + U^SaH - 2Uy^Hdt - 2U^a{d^ay - dyU^)) = 0, (C.3e) 


XB + (dA + V)XB j + Z^^4>bXb - ) If 


Maxwell r: 

Sd (^diSQ - - Z.^ado(p + /S {-f~^dodt + dxdx + dydy) + 2HUy^f^/‘^e~^‘^~‘^^^doT = 0, 

(C.3f) 


Maxwell i: 

dr (/-^e-2^d + S5* - + 2U^XBed'^d,ak = 

2U^aididi - di4)djT + f^l'^e^^^dAd^od - ^a") 

+ if/V2e('^-4)^a,- {6151 - 5l5i) - J:Si{5^Jl - d^Si) - ^Si{6l,6'; - S^S^a) ■ 

(C.3g) 


Note that 


denotes the totally antisymmetric symbol in flat space. 
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